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Abstract 



We present exact calculations of the Potts model partition function Z(G, q, v) 
for arbitrary q and temperature-like variable v on n-vertex strip graphs G of 
the triangular lattice for a variety of transverse widths equal to L vertices and 
for arbitrarily great length equal to m vertices, with free longitudinal boundary 
conditions and free and periodic transverse boundary conditions. These parti- 
tion functions have the form Z(G,q,v) = Y^j=f' X c z,G,j(^z,G,j) m_1 - We give 
general formulas for Nz : G,j an d its specialization to v = — 1 for arbitrary L. 
The free energy is calculated exactly for the infinite-length limit of the graphs, 
and the thermodynamics is discussed. It is shown how the internal energy cal- 
culated for the case of cylindrical boundary conditions is connected with critical 
quantities for the Potts model on the infinite triangular lattice. Considering 
the full generalization to arbitrary complex q and v, we determine the singular 
locus B, arising as the accumulation set of partition function zeros as m — > oo, 
in the q plane for fixed v and in the v plane for fixed q. Explicit results for par- 
tition functions are given in the text for L = 3 (free) and L = 3,4 (cylindrical), 
and plots of partition function zeros and their asymptotic accumulation sets 
are given for L up to 5. A new estimate for the phase transition temperature 
of the q = 3 Potts antiferromagnet on the 2D triangular lattice is given. 



Key Words: Potts model, triangular lattice, exact solutions, transfer matrix, Fortuin- 
Kasteleyn representation, Tutte polynomial. 
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1 Introduction 



The g-state Potts model has served as a valuable model for the study of phase 
transitions and critical phenomena PUI2|- m this paper we present some theorems on 
structural properties of Potts model partition functions on triangular-lattice strips of 
arbitrary width equal to L vertices and arbitrarily great length equal to m vertices. 
We also report exact calculations of Potts model partition functions for a number of 
triangular-lattice strips of various widths and arbitrarily great lengths. Using these 
results, we consider the limit of infinite length. For this limit we calculate thermo- 
dynamic functions and determine the loci in the complex q and temperature planes 
where the free energy is non-analytic. These loci arise as the continuous accumulation 
sets of partition-function zeros. This work is an extension to the triangular lattice of 
our earlier study for the square lattice [3]. 

Consider a graph G = (V, E), defined by its vertex set V and edge set E. Denote 
the number of vertices and edges as \V\ = n and \E\, respectively. For technical 
simplicity, we restrict to connected loopless graphs. On this graph G, at temperature 
T, the Potts model is defined by the partition function 

Z(G,q,v) = Y, e ~ m ( L1 ) 

with the (zero-field) Hamiltonian 

H = -jY,**<* i (1-2) 

(ij) 

where = 1, . . . , q are the spin variables on each vertex i 6 V; [3 — (fc^T) -1 ; and 
(ij) G E denotes pairs of adjacent vertices. We use the notation 

K = pJ, a = e K , v = a-l (1.3) 

so that the physical ranges are (i) a > 1, i.e., v > corresponding to oo > T > 
for the Potts ferromagnet, and (ii) < a < 1, i.e., —1 < v < 0, corresponding to 
< T < oo for the Potts antiferromagnet. One defines the (reduced) free energy per 
site / = —/3F, where F is the actual free energy, via 

f({G},q,v)= \im\n[Z(G,q,v)^ n } (1.4) 

where we use the symbol {G} to denote lim^oo G for a given family of graphs G. 

For our results in this paper we shall consider two types of boundary conditions: 
free and cylindrical. Here, free boundary conditions mean free in both the transverse 
and longitudinal directions (the latter being the one that is varied for a fixed width), 
while cylindrical boundary conditions mean periodic in the transverse direction and 
free in the longitudinal direction. Exact partition functions for arbitrary q and v have 
previously been presented for strips of the triangular lattice for width L = 2 in j3] 
and results have been given for L = 3, 4 in j^j. Here we shall give new explicit results 



3 



using the transfer matrix method (in the Fortuin-Kasteleyn representation (HI U\) 
for strips with free and cylindrical boundary conditions with width L = 3 (free) in 
section 3.2 and L = 3,4 (cylindrical) in sections 4.2 and 4.3. We have obtained new 
algebraic computer results for 4 < L < 6 (free) and 5 < L < 9 (cylindrical); these 
are too lengthy to include here but are included in the mathematica computer file 
transf er_Tutte_tri .m which is available with the electronic version of this paper in 
the cond-mat archive at http://arXiv.org We shall also present plots of partition 
function zeros and their asymptotic accumulation sets in the limit of infinite length, 
for widths 2 < L < 5. Where our results overlap with those for L = 2,3,4 given 
before in 0J|S], they are in agreement. As noted, we shall also give exact results valid 
for these strips with arbitrary width and length, namely the structural theorems of 
section 2. 

There are several motivations for this work. Clearly, new exact calculations of 
Potts model partition functions are of value in their own right. In addition, these 
calculations can give insight into the complex-temperature phase diagram of the two- 
dimensional (2D) Potts model on a particular lattice. This is useful, since the 2D Potts 
model has never been solved except in the q = 2 Ising case. From a mathematical 
point of view, the partition function of the Potts model on a graph G is equivalent to 
the Tutte polynomial on the same graph G (see below). Thus, we can extract very 
useful combinatorial information on the graph G. 

Let G' = (V, E') be a spanning subgraph of G, i.e. a subgraph having the same 
vertex set V and an edge set E' C E. Then Z(G,q,v) can be written as the sum 

Z(G, q ,v)=J2^ lE ' 1 (1-5) 

G'CG 

where k(G') denotes the number of connected components of G' . The formula (|1.5|) 
enables one to generalize q from Z + to R + (keeping v in its physical range), and it 
also shows that Z(G,q,v) is a polynomial in q and v (equivalently, a). 

The Potts model partition function on a graph G is essentially equivalent to the 
Tutte polynomial Ell HI] and Whitney rank polynomial [T21 H EEH1 El EES HUT 
Here the Tutte polynomial of an arbitrary graph G = (V, E) is 

T(G,x,y)=J2(x-±) k{G ' ) - kiG \y-±) ciG ' ) (1-6) 

G'CG 

where G' again denotes a spanning subgraph of G and c(G') denotes the number 
of independent circuits in G', satisfying c{G') = \E'\ + k{G') — \V\. Since we only 
consider connected graphs G, we have k(G) = 1. From (jl.5|) and it follows 

that the Potts model partition function Z(G, q, v) is related to the Tutte polynomial 
T(G, x, y) according to 

Z(G, q, v) = (x- l) kiG \y - 1)WT(G, x, y) (1.7) 

where 

x = 1 + - (1.8a) 

v 

y = a = v + l (1.8b) 
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so that 



q = (x- l)(y - 1) 



(1.9) 



In addition to the works previous exact calculations of Potts model partition 

functions for arbitrary q and v on lattice strips and/or studies of their properties 
include |3 [TTJ [TBI UHl 1201 123 I2U 123 a related early study of chromatic and 
Tutte polynomials for recursive families of graphs is [27] . 

Various special cases of the Potts model partition function are of interest. One 
special case is the zero-temperature limit of the Potts antiferromagnet, i.e., v = — 1. 
For sufficiently large q, on a given lattice or graph G, this exhibits nonzero ground 
state entropy So (without frustration). This is equivalent to a ground state degeneracy 
per site (vertex), W > 1, since So = ksliaW. The T = partition function of the 
g-state Potts antiferromagnet on G satisfies 

Z(G,q,-l) = P(G,q) (1.10) 

where P(G, q) is the chromatic polynomial (in q) expressing the number of ways of 
coloring the vertices of the graph G with q colors such that no two adjacent vertices 
have the same color [HI El EEl I2H| • The minimum number of colors necessary for this 
coloring is the chromatic number of G, denoted x(G). We have 

W{{G},q) = lim P{G,q) 1/n (1.11) 

In the context of our current work we recall that the chromatic number for the 2D 
triangular lattice is x(tri) = 3. This chromatic number also applies to strips of 
the triangular lattice with free longitudinal boundary conditions and free transverse 
boundary conditions. For the triangular-lattice strips with cylindrical (i.e., free longi- 
tudinal and periodic transverse) boundary conditions, x = 3 if the width L = mod 
3 and x = 4ifL = lor2 mod 3. References to papers on the special case v = — 1 
are given, e.g., in [TBI E01 E] • 

Using the formula (jl.5|) for Z(G,q,v), one can generalize q from Z + not just to 
R + but to C and v from its physical ferromagnetic and antiferromagnetic ranges 
< v < oo and — 1 < i; < to w 6 C. A subset of the zeros of Z in the two-complex 
dimensional space C 2 defined by the pair of variables (q, v) can form an accumulation 
set in the n — > oo limit, denoted £>, which is the continuous locus of points where 
the free energy is nonanalytic. This locus is determined as the solution to a certain 
{G}-dependent equation. For a given value of v, one can consider this locus in the q 
plane, and we denote it as B q ({G}, v). In the special case v = — 1 where the partition 
function is equal to the chromatic polynomial, the zeros in q are the chromatic zeros, 
and B q ({G},v = —1) is their continuous accumulation set in the n — > oo limit. With 
the exact Potts partition function for arbitrary temperature, one can study B q for 
v 7^ — 1 and, for a given value of q, one can study the continuous accumulation set 
of the zeros of Z(G, q, v) in the v plane (complex-temperature or Fisher zeros [STJ - 
other early references include [321 EHl Ell ) ■ This set will be denoted B v ({G},q). 
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2 General Results for Recursive Families of Graphs 



A recursive family of graphs is one in which one constructs successive members of 
the family in a recursive manner starting from an initial member. Recursive families 
of graphs that are of particular interest here are strips of regular lattices of a given 
width L vertices (with free or cylindrical boundary conditions) and arbitrarily great 
length m vertices (with free boundary conditions). 

A general form for the Potts model partition function for the strip graphs consid- 
ered here is [TH] 

Z(G,q,v)= ^,(A,, ; r ' (2-1) 

3=1 

where the coefficients cqj and corresponding terms Xqj, as well as the total number 
Nz,g,\ °f these terms, depend on the type of recursive graph G (width and boundary 
conditions) but not on its length. (In ^H], a slightly different labelling convention 
was used so that Xqj rather than A^J 1 appeared in the summand of eq. (|2.ip .) In the 
special case v — — 1 where Z reduces to the chromatic polynomial (zero-temperature 
Potts antiferromagnet), eq. (j2.1|) reduces to the form [HH] 

P(G,q)= cajiXp^r- 1 (2.2) 

3=1 

For the lattice strips of interest here, we define the following explicit notation. Let 
-^z,sq,Bc t Bc e ,L,\ denote the total number of A's for the square-lattice strip with the 
transverse (t) and longitudinal (£) boundary conditions BC 4 and BC^ of width L. 
Henceforth where no confusion will result, we shall suppress the A subscript. The 
explicit labels are N z , S q,FF,L and N z ,tri,FF,L for the strips of the square and triangular 
lattices with free boundary conditions, and N z , S q,PF,L and N z ^ti,vf,l for the strips of 
these respective lattices with cylindrical boundary conditions. 

For the lattice strip graphs of interest here we can express the partition function 
via a transfer matrix T (in the Fortuin-Kasteleyn representation) of fixed size M x M: 

Z(G,q,v)=tr[A(q,v)-T(q,v) m - 1 ] (2.3) 

which then yields the form ()2.1|) . Since the transfer matrix T and the boundary- 
condition matrix A are polynomials in q and v, it follows that the eigenvalues {A^} 
of the transfer matrix and the coefficients cqj are algebraic functions of q and v. 

One of the basic structural properties of the Potts model partition function on a 
given strip is the number of different eigenvalues of the transfer matrix (in the Fortuin- 
Kasteleyn representation), Nz,g,bc,l, in eq. f|2.1|) . In in addition to proving various 
formulas for these numbers for certain strip graphs, we presented a conjecture (de- 
noted Conjecture 3 in j3] and given as Theorem 4.3.5 by one of us (S.-C.C.) in [36j). 
We now give the proof of this result. 
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Theorem 2.1 For arbitrary L, 




iV Z,sqtri,FP,| 
1 AT 

2 iV Z,sqtri,FP,^ 



for even L 
for odd L 



(2.4) 



-/Vz,sqtri,FP,L 



2L 



L 



) 



(2.5) 



Proof We recall first that the quantity 2Nz )S q,FF,L — Nz,tn,FF,L discussed in |3] 
gives the number of non-crossing partitions for a transverse slice of these two respec- 
tive strips (which is the path graph Tl) such that these partitions are symmetric under 
reflection about the longitudinal axis. The quantity 2Nz, S q,PF,L — ^,tri,PF,L gives the 
corresponding number of non-crossing partitions for a transverse slice (which is the 
circuit graph Cl) of the two respective cylindrical strips such that these partitions 
are symmetric under reflections about the longitudinal axis and rotations around this 
axis (the latter being included since there is no special azimuthal direction). We shall 
prove eq. (|2.4j) for odd L first and then for even L. 

For odd L, let n = (L+l)/2, and denote 2iV ZiSqi PF,L=2ri-i-A r z,tri,PF,L=2n-i as X n for 
simplicity. Consider a transverse slice with periodic boundary conditions. Since this 
is topologically invariant under rotations around the longitudinal axis, we can label 
one vertex as 1 and other vertices 2,3, ... ,n, ml , (n — 1)', . . . , 2', in a counterclockwise 
manner (relative to a specified longitudinal direction), and consider the reflection 
symmetry with respect to the longitudinal axis passing through vertex 1. In order to 
classify the types of colorings of the vertices, we shall introduce diagrams consisting 
of the L vertices on this transverse slice. In this context, we shall refer to two vertices 
as being "connected" if these have the same color and shall denote this by using the 
Kronecker delta function. 

The sets Px n of partitions for n = 1, 2, 3 that are invariant under this reflection 
symmetry are: P Xl = {!}, Px 2 = {1, £2,2', ^1,2^1,2' = £1,2,2'}, and 



We can classify these partitions into cases that have m vertices on one side of the 
slice (including vertex 1) connected to at least one other vertex (on the same side 
or the other side, indicated by the primes above), where < m < n. For m — 0, 
this is the partition 1, that is the identity partition, defined as the one in which 
all blocks are "singletons", i.e., there are no connections in the sense given above. 
For m = l, there is only one possibility: S xltX ' , where 2 < x\ < n. For m = 2, 
let us denote the connected vertices clS X2 and, with no loss of generality, take 
X\ < X2, then there are the partitions ^ x ^,x 2 -. &x 1 ,x 2 ,x' 1 i and 5 x 1 ,x' 1 ^x 2 ,x' 2 - F° r t ne first 
case, x\ can be vertex 1, but for the second and third cases, 2 < x\ < n since x\ 
and x\ are different vertices. The corresponding partitions 8 X > iX / for the first and 
second cases can be obtained by reflection symmetry from the ones that we have 



{1,(52,2', £3,3', £l,2,2', £l,3,3', £2,3£2',3', £2,2'£3,3' 
£2,3,2',3', £l,2,2'£3,3', £l, 2,3,2', 3' } 



(2.6) 
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listed and hence, for simplicity, are not shown. With no loss of generality, we take 
xi < x 2 < . . . < x m . For m = 3, the partitions are <5 X1 , X A 3 ,4> 5 X1 , X2 , X3 , tf^^y^y, 
S X1 , X2 , X3 , X ' V tfxi.zi&w.xsj ^x u x'^x 2 ,x 3 ,x' 2 , ^1,^^2,4^3,4- Having given these illustrations 
of the specific partitions for < m < 3, we next proceed to the general case. 

The partitions that have m vertices on one side of the slice connected to as least 
one other vertex (on the same or opposite side) can be classified further. Let us 
denote a m as the number of cases where the vertex x\ has only the connection S xliX i 
for 1 < m. In these cases, x\ cannot be vertex 1, and the number of the partitions 
for each m < n — 1 is ( n_1 ) • The last three partitions of m = 3 given above are 
examples of these cases. We denote b m as the number of cases where the vertex x\ 
has connection to at least one unprimed vertex with or without 5 xi y for 2 < m. The 
first and second partitions for m = 3 are examples of these cases without 6 XliX > , and 
the third and fourth partitions for m = 3 are examples of these cases with S Xl>x i . 
Notice that while X\ cannot be vertex 1 for the cases with 5 XljX ' , which have the 
number of partitions ( n ~ 1 ) for each m < n — 1, xi can be vertex 1 for the cases 
without 5 Xlt x i , and the number of partitions for each m < n is (^). Therefore, 

m=l V 7 m=2 v 7 m=2 v 7 

Next, we shall obtain expressions for a m and b m . The cases for be obtained 

from all the cases with m — 1 vertices by changing X{ to Xi + i for 1 < % < m — 1 and 
adding 5 xljX ' , i.e., a m is the same as the total number of the cases with m — 1 vertices, 

a m = a m _i + 26 m _i (2.8) 

Now consider the cases for b m without 5 Xl y . These can be further divided into two 
possibilities: the cases with 5 XltX2 and the cases with S XljXi where 2 < i < m. Denote 
the numbers of these two possibilities as d m and e m , respectively. Clearly, 

bm = d m + e m (2-9) 

The cases for d m can be obtained from the cases for 6 m _i by changing Xi to Xi+i for 

1 < % < m — 1 and adding 5 XUX2 , and from all the cases with m — 1 vertices, where 
the number is a m _i, by changing Xi to x i+2 for 1 < i < m — 2 and adding S xljX2 . 
Therefore, 

d m = a m _i + 6 m _i (2.10) 

The cases for e m can be obtained from the cases for rf m _ 2 by changing Xi to for 

2 < « < m — 2 and adding 5 X2iX3 , from the cases for c? m _ 3 by changing to Xi + % for 
2 < i < m — 3 and adding S X2jX3jX4 , etc. In general, the cases for e m are obtained 
from the cases for d m ^i by changing Xi to Xi + e for 2 < i < m — £ and 2 < £ < m — 2, 
and adding a set of £ vertices, where all the vertices in this set must have at least 
one connection to at least one other vertex in this set without reflection or rotation 
symmetry. The number of all possible connections of this set with £ vertices is the 
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Riordan number (given as sequence A005043 in |H7|), which will be denoted as rg and 
may be defined via the generating function |3£ 



R(z) 



l + z -( 1 _2z-3z 2 ) 1 / 2 ^ 

— = 2^ riZ 



2z(l + z) 



Therefore, 



m-2 
1=2 



-i r e 



In terms of generating functions 

A{z 



B(z 
D(z 
E(z 
R(z 



d\Z + a 2 z 2 + a 3 z 3 + 



b 2 z 2 + b 3 z 3 + b 4 z 4 + ... 
(I2Z 2 + d 3 z 3 + d^z 4 + . . . 



e 4 z 4 + e 5 z 5 + e 6 z 6 + . . . 



1 + r 2 z + r 3 z 3 + r 4 z + 



we can re-express eqs. (j2.8|) to (j2.12|) as 

A(z) — a\z 

z 

B{z) = 
D(z) 



z 

E(z) 



A(z)+2B(z) 
D{z) + E{z) 

A{z)+B{z) 
D(z)(R(z)-l) 



We find 



A(z) 



VI -2z-3z 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



The coefficients in the expansion of A(z) in eq. f)2. 13j) are, up to a shift, the central 
trinomial coefficients (given as sequence A002426 in [27]), i-e., for each value 1 < m 
the largest coefficient of (1 + z + z 2 ) m_1 . Next, 



B(z) 



1 - z 



2 l v /l-2z-3z 2 



- 1 



(2.16) 



The coefficients in the expansion of B(z) in ()2.13j) are given by the coefficients for 
the next-to-central column in the expansion of (1 + z + z 2 ) m_1 for 1 < m (listed as 
sequence A005717 in [HTJ). For D(z) we have the closed form 



D(z) 



l + z 



2 LV 1 -3z 
9 



(2.17) 



The coefficients in the expansion of D(z) in ()2.13|) are the numbers of directed animals 
of size m — 1 on the square lattice for 1 < m (given as sequence A005773 in |37j). 
Finally, we have 



E(z) = z 



1 - 3z + (1 - 2z)y/l-2z-3z 2 
l-2z-3z 2 + (3z - l)y/l-2z-3i 



l + z 



(2.18) 



The coefficients in the expansion of E(z) in ()2.13j) are the coefficients forming the 
second column from the center, in a tabular format, in the expansion of (1 + z + z 2 ) m ~ 2 
for 2 < m (given as sequence A014531 in |37j). 

Recall the binomial transformation for two sequences of numbers [sq, s\, s 2 , ■ ■ ■], 
[to,ti,t 2 , ■ ■ ■] with generating functions S(z) = Y^=o s nZ n and T(z) = J2'^=o t nZ n - If 
these sequences have the relation 



n 

\ / n 



m=0 



m 



then 



T(z) 



1 



In our case, if the generating function of X n is X{z) 
combine eqs. I|27j) . (j2~TH|) and (gHED to get 



(2.19) 
(2.20) 

Y^=i X n z n , then we can 



X(z) 



z 



1 - z 



VT^Iz 



1 - 
- 1 



-.4 



1-z 



+ 



-B 



1-z 



+ 



1 - z 



-B 



This has the expansion 



Xtz) 



E 

n=l 



2n - 1 



n 



(2.21) 



(2.22) 



Note that ( 2n ~ l ) = i( 2 ™) = i( (L +t)/ 2 ) (g iven as sequence A001700 in jSTj). 

Let us next proceed to consider the number of non-crossing partitions for a slice 
of the transverse vertices which has periodic boundary conditions and reflection sym- 
metry for even L. Denote n = L/2, and 2N z , S q,PF,L=2n — ^z,tri,PF,L=2n by Y n for 
simplicity. There are two possibilities: the reflection axis does not go through any 
vertex or goes through two vertices. These possibilities will be denoted as type I 
and and type II partitions, respectively, and the number of partitions of each of 
these two types as and Y* 1 . For type I partitions, label the vertices on one 
side of the reflection axis as 1, 2,. . . , n and the corresponding reflected vertices as 
1', 2', . . . , n' on the other side. For type II partitions, label the vertices as 1, 2, 
. . . , n, n + 1, n', . . . , 2', where vertices 1 and n + 1 are on the reflection axis. 
The sets Py^i of type I partitions for n = 1,2 having this reflection symmetry are: 
P y i = {1,(51,1'} and V Y i = {1, S^y, 8 2 ,2', hpSvp', ^1,2,1', t^i/f^'}- The sets Pyn of 
type II partitions for n = 1,2 having this reflection symmetry are: P y n = {1,(^2} 
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and Pyii = {1, $2,2', £l,3, ^1,2,2', ^2,3,2', ^i,2,3,2',3'}- We notice that partitions 1, i.e. iden- 
tity and Si : 2,...,l (i-e., a unique block) are both contained in the type I and type II 
classes of partitions. Since we have rotational symmetry, the partitions which are not 
symmetric with respect to the central axis perpendicular to the reflection axis are 
counted twice in either type I or type II classes of partitions. A similar statement 
applies to the partitions which are symmetric with respect to the axis perpendicular 
to the reflection axis if L is a multiple of 4. If L is not a multiple of 4, the partitions 
that are symmetric with respect to the perpendicular axis are counted once in both 
type I and type II classes. Therefore, 2Y n is the sum of all possible partitions in these 
two classes of partitions. 

We again classify these partitions into cases which have m vertices on one side 
of the slice connected to at least one other vertex (including vertex 1 and vertex 
n + 1 for type II partitions). It is clear that < m < n for type I partitions and 

< m < n + 1 for type II partitions. Let us consider type I first. For m = 0, there 
is the partition 1 (identity). For m — 1, there is only one possibility: 5 X1)X ' , where 

1 < x\ < n. For m — 2, there are S XljX2 , 5 xl)X2iX ' and S X1;X '6 X2;X ', where x\ < x^. 
The corresponding 5 X > )X ' for the first and second cases are not shown for simplicity. 
We shall again take X\ < x<i < . . . < x m . The number of the cases with m vertices 
connected is a m + 2b m = a m+1 as discussed above eq. (|2.8|) . and the number of the 
partitions is ( n ) for each < m < n. Therefore, 



Let the generating function of Y* be Y 1 (z) = Y^=i^n zn - Using eq. (|2.15|) and 
modifying eq. (|2.2(J|) . we have 



For type II partitions, neither vertex 1 nor vertex n + 1 has a corresponding 
symmetric partner, 1' and (n + 1)', respectively. We have to classify the cases where 
Xi is only connected to at least one other unprimed vertex (in these cases, x\ can be 
vertex 1), the number of which cases was denoted b m earlier, into two possibilities: the 
cases where x m is only connected to at least one other unprimed vertex (so that x m 
can be vertex n + 1), and the cases where x m is connected to x' m (so that x m cannot 
be vertex n + 1) among other possible connections. The number of partitions for these 
two possibilities will be denoted as f m for 2 < m and h m for 3 < m, respectively. The 
partitions where x\ is connected to x[, i.e., 5 XlyX ' v among other possible connections 
can also be classified into two possibilities: the cases where x m is only connected to 
other unprimed vertices, and the cases where x m is connected to x' m among other 
possible connections. The number of partitions for these two possibilities will be 
denoted as h' m for 3 < m and i m for 1 < m, respectively. Notice that h' m = h m 
because of the reflection symmetry. We shall need another set of partitions where X\ 
is connected to x m with possible connection to other unprimed vertices, and all the 




for 1 < n 



(2.23) 




(2.24) 
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other vertices are connected to at least one other unprimed vertex; we denote the 
number of partitions for them as j m for 2 < m. If we add an additional connection 
&n,a;i to these partitions, and denote the number of partitions as k m , then k\ = 1 and 
& m = jm for 2 < m. 

For the cases for f m , x\ and x m can be vertices 1 and n + 1, respectively, so the 
number of the partitions is ("j^ 1 ) for each m. For the cases for h m , x% can be vertex 
1 but x m cannot be vertex n + 1, so the number of the partitions is ( n ) for each m. 
For the cases for i m , neither x\ can be vertex 1 nor x m can be vertex n + 1, so the 
number of the partitions is ( n ~ 1 ) for each m. Therefore, 

f^+E^'l^^t)^^ 1 ) (2 ' 25) 

m=2 v 7 m=3 v 7 m=l v 7 

We next obtain expressions for f m , h m and i m . From the definitions, 

f m + h m = b m (2.26) 

and 

h' m + i m = a m + frm = ^m+i j (2.27) 

where we use eq. (|2.10p . The cases for f m include all the cases for j m , plus the cases 
for both f m -2 and h m -2 with the additional connection S Xm _ ljXm , plus the cases for 
both f m -3 and h m _ 3 with the additional connection S Xm _ 2tXm _ uXm , etc. We have 

m—2 m— 2 

fm = jm + 2_j 3t lfm-e + h m -l] = j m + j<? &m-£ (2.28) 

€=2 1=2 

The cases for fo m can be obtained by adding S XmtX i to the cases for f m -i and /i m _i, 
and adding S Xm _ 1>Xm>x r m to the cases for j m -2 and /i m _2, etc. We have 

m— 2 m— 2 

— ^3 ^ [/m-^ + h-m-l] = ^3 ^ &m-^ (2.29) 
1=1 1=1 

The cases for h' m can be obtained by adding 6 Xm _ 1)Xm to the cases for h' m _ 2 and z m _2, 
and adding S Xm _ 2)Xm _ ljXm to the cases for h' m _ 3 and i m _3, etc. We have 

m— 1 m— 1 

= 2 J/ [^m-/ + imr-t] = ^Jt d m -i+l , (2.30) 
1=2 1=2 

which should be equal to h m , as mentioned before. The cases for i m include all the 
cases for k m , plus the cases for both hl m _ x and % m -\ with additional connection S Xm)X i , 
plus the cases for both h' m _ 2 and i m -2 with additional connection S Xm _ uXmjX r m , etc. We 
have 

m— 1 m— 1 

= k m + fcf + i m -i\ = km+^kg d m -i+i (2.31) 

i=i i=\ 
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The cases for j m can be obtained from the cases for r m by removing the cases for r m _2 
with additional S Xm _ 1)Xm , removing the cases for r m _ 3 with additional S Xm _ 2tXm _ 1>Xm , 
etc. We have 

m-2 

jm = r m-^2 I* Tm - e {222) 
1=1 

In terms of the generating functions 

F(z) = f 2 z 2 + f 3 z 3 + hz 4 + ... 

H(z) = H{z)' = h 3 z 3 + h 4 z 4 + h 5 z 5 + . . . 

I(z) = i 1 z + i 2 z 2 + i 3 z 3 + . . . 

J{z) = j 2 z 2 + j 3 z 3 + ]A z A + . . . 

K(z) = hz + k 2 z 2 + k 3 z 3 + . . . (2.33) 
we can reexpress eqs. (I2.26j) to (|2.32|) as 

F{z) + H(z) = B{z) 

H(z) + I(z) = ^ 

z 

F{z) = J{z) + B{z)J(z) 

H(z) = B{z)K{z) = ^^J{z) 

z 

I(z) = K{z) + ^^K{z) 



We find 



J(z) = R{z)-l-(R(z)-l)J(z) = (R(z)-l){\- J{z)) (2.34) 



R(z) 



J{z) = = ^-z- VI - 2z - ^) , (2.35) 

which is essentially the Motzkin number [101 E] (given as sequence A001006 in [37J). 
This can be understood since the cases for j m are in one-to-one correspondence with 
the non-crossing, non-nearest-neighbor partitions of m— 1 vertices with free boundary 
condition, as given in eq. (2.1.1) of [3], i.e., j m = M m _ 2 . We also have 

K(z) = z + J{z) = -(1 + z - Vl - 2z - 3z 2 ) (2.36) 
The first five equations in eq. 1)2.34)1 are redundant, and can be solved to have 

F{z) = f = zA(z) , (2.37) 

V 1 — 2z — 3z z 
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i.e., f m = a m -i, which is the largest coefficient of (1 + z + z 2 ) m 2 for 2 < m. Next, 
we have 

'1 - 2z 



Hz) 



1-32 



\A-2z-3^ 2 



1 



(2.38) 



The coefficients in the expansion of H(z) are the numbers of directed animals of size 
m — 2 on the square lattice with the first two quadrants (which can grow in right, left 
and up directions) for 2 < m (given as sequence A005774 in pTj). Finally, 



/(:) = : J ±±±- = 2D(z) + z 



(2.39) 



which is given as sequence A0025565 in [37]. 

Denote the generating function of Y* 1 as Y 11 (z) = Y^=i^n zn - Combining eqs. ([2.25)1 . 
1)2.37}) . (|2.38|) . ()2.39|1 . and using the binomial transformation ()2.20|) . we obtain the re- 
lation 



Y n (z) 



1 



l-z 
1 

l-z 



1 — 7T~~ ) + T— H{ 



z(l — z) 1 — z 1 — z 1 — z 1 — z 1 — z 



z 



+ 



1 



1 -3z 



1-2) 2 V / T^4i 2(l-^) 2 L v /r^4^ 

2 



+ z-l 



z 



[l-zf^l^Ai 



1 



- 1 



which is the same as the Y l (z) in eq. (|2.24|) . We finally have 

Y l (z) + Y u (z) 1 



(2.40) 



Y z) 



JT=I~z 



- 1 



(2.41) 



which is the generating function of ( ™) = ( L / 2 ), i.e., the central binomial coefficients 
given as sequence A000984 in |37j . 

Eqs. fl2.21j) and (|2.41|l prove the theorem for odd and even L, respectively. □ 

Another structural result was given as Conjecture 4 in |S] and Theorem 4.3.7 in 
[36 for the number of A's for the triangular lattice with cylindrical boundary condition 
for the Potts model partition function. We state this and give the proof. 



Theorem 2.2 For arbitrary L, 



N 



Z,tri,PF,L 



c L + Yl ^ L / d ) 

d\L; l<d<L 



2d 

d 



(2.42) 



where d\L means that d divides L and 4>(n) is the Euler function, equal to the number 
of positive integers not exceeding the positive integer n and relatively prime to n. 
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Proof As shown in , Nz,tii,FF,L = Cl is the number of non-crossing partitions 
of a set of L vertices with free boundary conditions. LNz t tn,PF,L ~ ^Yz,tri,FF,z is equal 
to 2(L — 1) for prime L, as shown in [HUSH], since all the partitions of A^ j t r i,pF,£ have 
periodicity L for prime L except for the partitions 1 (i.e., all blocks being singletons) 
and 5i,2,...,l (i-e., a unique block) which have periodicity 1. Consider a general L and 
assume that L has the factor d (so that for prime L, d can only be 1 or L). Denote 
the number of partitions which have periodicity d modulo rotations as 2o^. Then 

LN ZM , PFtL - iVz, t ri,FF,L = Yl 2a ^ L ~ d ) ( 2 ' 43 ) 

d\L 

where d are all the positive integers that divide L. 

There are two kinds of partitions among 2oid, each of which contains specific 
partitions. The first kind of partitions is defined by the condition that each set of d 
adjacent vertices does not have connection with any other vertex in the complemen- 
tary subset of vertices. Let us label the vertices in one set of d vertices as 1,2,. . . ,d. 
For d = 1, the set of partition is {1}. For d = 2, the set of partition is {5i >2 }. For 
d = 3, the set of partitions is {£1,2, £1,3, ^1,2,3}- The second kind of partitions is de- 
fined by the condition that all of the first vertices of each of the sets of d vertices are 
connected to each other. For d — 1, the set of partition is {£1,1'}, where 1' is the first 
vertex of the adjacent set of d vertices. For d = 2, the set of partition is {#1,1'}. For 
d = 3, the set of partitions is 5i,i'#2,3, ^i,i',3}- Since there is a one-to-one corre- 

spondence between these two kinds of partitions, let us only consider the second kind 
of the partitions which always have the connection S^y. Denote a' d as the number of 
partitions which has periodicity less than or equal to d modulo rotations. This can 
be written in term of as 

ai = 5>* (2-44) 

d'\d 

For the set of d vertices, if vertex 1 does not connect to any other vertex except V, the 
first vertex of an adjacent set, then the number of partitions is equal to iYz,tri,FF,d-i — 
Cd-\- If vertex 1 has a connection to vertex b for 1 < b < d in addition to vertex 
1', then the number of partitions is C^Cd-b- The connection of vertex 1 to vertex 
d—b+2 is equivalent to the connection of vertex 1 to vertex b under rotation symmetry, 
and should not be considered again. Of course, vertex 1 can have connection to more 
than one vertices in addition to 5i^'. Therefore, to calculate a' d , we have to partition d 
first. For each partition of d, e.g., [xi,x 2 , ■ ■ ■ , xj\ with x\ + x 2 + ■ ■ ■ +Xi = d, where i is 
the numbers of vertices in the set of d vertices connected to vertex 1 (including vertex 
1), we multiply all the corresponding C Xi -iS by the number of different combinations 
of this partition modulo rotations (thus [2,2,1,1] is the same as [2,1,1,2] and [1,1,2,2], 
but different from [2,1,2,1]). As an example, for d = 5, we have the partitions 
5 = 4 + 1 = 3 + 2 = 3 + 1 + 1 = 2 + 2 + 1 = 2 + 1 + 1 + 1 = 1 + 1 + 1 + 1 + 1, and 
therefore 
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= 14 + 5xl + 2xl + 2xl 2 + l 2 xl + lxl 3 + l 5 = 26 (2.45) 

The factors |, |, and \ are included because of the equivalence under rotations. 

We find that a' d is the same as the rooted planar trees with d edges with the 
property that rotations about the root vertex yield equivalent trees, given as sequence 
A003239 in |37|. The reason can be explained as follows. Define a "planted" tree as a 
rooted tree with the property that the root vertex has degree one (where the degree 
of a vertex is the number of edges connected to it). Now we can construct planted 
trees from subtrees of the rooted planar tree that are each connected to the root 
vertex. Note that the number of planted tree with d edges is Cd-i [SB]- The sequence 
A003239 also gives the number of necklaces with a total of 2d beads where d beads 
have one color and the other d beads have another color. The explicit formula is 



Al I J \ ' / 



a 

•2d 

d'\d 

where <f>(n) was defined above after eq. (|2.42|) . 

Now the number of partitions that have periodicity equal to d modulo rotations, 
aa, is the Mobius transformation [39J of a' d , given by 

a d = ^n{d/d')a d , (2.47) 

d'\d 

where //(n) is the Mobius function, defined as —1 if n is prime, if n has a square 
factor, and 1 for other n. We find that ad is listed as sequence A022553 in [37|, which 
is the Lyndon words containing a total of 2d letters with d letters of one type, and 
the other d letters of another. One has the explicit formula 



d'\d ^ ' 

i.e., 2dad, is the Mobius transformation of [ 2 d ), or equivalently, 

' d'\d 

Therefore, the total number of partitions that have periodicity d modulo rotations is 

2«* = 3 £>( d A0 ( d?) (2-50) 

d'\d ^ ' 

This is given as sequence A060165 in [3*7j . 
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We finally have 
LN Z 

,tri.PF,L -N z .tri.FF.L 



^2a d (L-d) 

d\L 

2^ 2La. 

d\L 



e^e^')( 2 ;)-( 

rf|L d'\d v 7 v 

E E 

d'\L d\L;d'\d 



2L 



2L 



5>(L/<0 



d'\L 



2d! 



2L 

d' ) \ L 



d\L; l<d<L 



d 



(2.51) 



and the theorem follows. □ 



Theorems 12.11 and 12.21 imply a simple corollary which gives an exact formula for 

Nz,sq,PF,L- 



Corollary 2.1 For arbitrary L, 

Nz,sq,PF,L 



Proof For even L is trivial. For odd L, we only have to notice that 
L + l \ L(L — l)--f _ f L \ ( L \ 

~~ (L - i)A 



(1/2) 



XL + l)/2 
This proves the corollary. □ 



' L-\ - 
. 2 - 



[L/2}J 



(2.52) 



(2.53) 



One of us (S.-C.C.) has presented a conjecture, Conjecture 4.2.2 in j^Hj, which 
generalizes Theorem 3 in |3j (i.e., Theorem 4.2.6 in [HHl)- We restate this conjecture 
here. Consider the number of A's in the chromatic polynomial for a strip of the 
triangular lattice with cylindrical boundary conditions: 



Conjecture 2.1 For arbitrary L, 



N. 



1 



P,tri,PF,L 



(2.54) 

d\L; l<d<L 

where td = dd+i is the largest coefficient in the expansion of (1 + z + z 2 ) d , i.e., the 
central trinomial coefficient (given as sequence A002426 in \3 r lj ), and di is essentially 
the Riordan number r^ (given as sequence A005043 in 

1 L = 1 

r L L>2 



di 



(2.55) 
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where the reader is cautioned not to confuse the di in eq. (|2.54|) with the different 
quantity d m in eqs. (|2.9|) and (j2.10j) . We can motivate this conjecture as follows. 
We know that di is the number of non-crossing non- nearest-neighbor partitions of a 
set of L vertices with periodic boundary conditions, as shown in [T^]. The number 
LNp t tri,PF,L — di is equal to L — 1 for prime L, as shown in [3J EE], since all the 
partitions of iVp.tri.PF.z, have periodicity L for prime L except for the partitions 1 (i.e., 
all blocks being singletons) which have periodicity 1. Consider a general L having 
the factor d, and denote the number of non-nearest-neighbor partitions that have 
periodicity d modulo rotations as (3d- Then 

LN PM ^ L -d L = Y,Pd(L-d) (2.56) 

d\L 

There are two kinds of partitions among (3d] we denote the number of these as (3\ 
and (3 d l , respectively. The first kind of partitions is defined by the condition that each 
set of d adjacent vertices does not have connection to any other vertices in other sets. 
Label the d vertices as 1,2,. . . ,d. For d = 1, the set of partitions consists of just {1}. 
There is no non-nearest-neighbor partition for d = 2. For d = 3, the set of partitions 
is {^1,3}- For d = 4, the set of partitions is {61^,61^}. The second kind of partitions 
is defined by the condition that all of the first vertices of each set of d vertices are 
connected to each other. There is no non-nearest-neighbor partition for d — 1. For 
d = 2, the set of partitions is comprised of {£1,1'}, where 1' is the first vertex of an 
adjacent set of d vertices. For d = 3, the set of partitions is For d = 4, the set 

of partitions is {£1,1/, <5i,i'<52,4}- Notice that for d = 4, {<5i,i',3} is the same as {£1,1' } 
for d = 2, and should not be included to avoid double-counting. 

Consider the first kind of partitions. If the first vertex 1 is connected to the 
last vertex d among other possible connections with other vertices in the set of d 
vertices, the number of these partitions is rd-i, the Riordan number (given as sequence 
A005043 in [HI])- This can be seen by identifying the vertices 1 and d to form a circuit 
with d — 1 vertices [19j . In addition, the set of d vertices can also be partitioned into 
several parts where the first and the last vertices of each part are connected, but 
parts are not connected to each other. To calculate (3 d , we partition d, and apply 
the same transformation as illustrated in Theorem 12.21 on r^-i, then apply a Mobius 
transformation. The first few values of (3 l d are 1, 0, 1, 2, 5, 11, 28, 68, 174, 445 for 
1 < d < 10. 

Consider the second kind of partitions. By an argument similar to that given in 
the proof of Theorem 12.21 if vertex 1 is not connected to any other vertex except V 
among the set of d vertices, then the number of partitions is just iVp^n^F^-i = Md-2- 
In addition, the vertex 1 can also be connected to other vertices among the set of d 
vertices. To calculate (3 d l , we again partition d, and apply the same transformation, 
as illustrated in Theorem 12 .21 on Md-2- We then apply a Mobius transformation. The 
first few values of 0% are 0, 1, 1, 2, 5, 11, 28, 68, 174, 445 for 1 < d < 10. We find 
that (3\ an d Pd are the same up to d = 10, except for d = 1 and 2. Although (3d is not 
listed in [37j, we find that the following relation holds up to d = 10, and is similar to 
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eq. JIHID : 

P d = \^n{d/d!)t d (2.57) 



d'\d 



If this is correct for arbitrary d, then the conjecture is proved. 
Combining Conjecture 12. II and Conjecture 1 in jHj, we have 



Conjecture 2.2 For arbitrary L, 

: N.P,tri,PF,£ + |^P, sq tri,FP,f 
: A^ Pitri ,PF,L + I^P, sq tri,FP,i±i ~ \ T ^ 



-^P,sq,PF,L 



/or ei>en L 
for odd L > 3 



(2.58) 



where iVp jSq tri,FP,i was given zn and r n z's i/ie Riordan number (given as sequence 
A0Q5Q43'in§]j). 

3 Potts Model Partition Functions for Triangular- 
lattice Strips with Free Boundary Conditions 

The Potts model partition function Z(G, q, v) for a triangular-lattice strip of width 
L and length m with free boundary conditions is given by 

Z(L F xm F ,q,v) = v T • H • T m ~ l ■ u id (3.1a) 
= w T -T m - 1 -u id (3.1b) 

where w T = v T ■ H. Hereafter we shall follow the notation and the computational 
methods developed in [3 Ell H21 HH] (For chromatic polynomials, a related matrix 
formulation has been discussed in j^HU|^l|12]-) Concerning notation, no confusion 
should result between the vertex set V, the variable v, and the vector v. Here 
T = V ■ H is the transfer matrix, and H (resp. V) corresponds to adding one more 
layer of horizontal (resp. vertical) bonds. The matrices T, V, and H act on the 
space connectivities of sites on the top layer, whose basis elements are indexed v-p 
are indexed by partitions V of the single-layer vertex set {1, . . . , L}. In particular, 

u id = V {{ i }i{2 },...,{L}}. 

To simplify the notation, we shall denote the elements of the basis v-p by a short- 
hand using Kronecker delta functions: for instance v/ji 31 /2},{4,5}} wm be written 
^1,3^4,5- We denote the set of basis elements for a given strip as P = {v^}. 

As we are dealing with planar lattices, only non-crossing partitions occur. The 
number of such partitions is given by the Catalan numbers 

C n — — 7— ( ^ \ (3.2) 



n + 1 V n 
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In the triangular-lattice strip with free boundary conditions there is no additional 
symmetry that allows us to reduce the number of partitions; thus the dimension of 
the transfer matrix for width L is 

N z ,tri,FF,L = C L (3.3) 

We have obtained the transfer matrices T(L-p) and the vectors w and u ic j using 
symbolic computation with mathematica as in 0]- We have double-checked 
these results using a different program written in perl. 

An equivalent way to present a general formula for the partition function is via a 
generating function. Labelling a lattice strip of a given type and width as G m , with 
m the length, one has 

oo 

T(G,q,v,z) = ^2z m Z(G m ,q,v) (3.4) 



with 



m=0 

where T(G, q, v, z) is a rational function 

TVr» ^ N(G,q,v,z) 

r(G, q, v, z) = — — (3.5 

V{G,q,v,z) 

dc gz (M) 

Af(G,q,v,z) = J2 A GJ Z ' ( 3 - 6 ) 

3=0 

Nz,tri,BC,L N Z} tri,BC,L 

V(G,q,v,z) = l+ J2 h G,3 z3 = II C 1 -^) ( 3 - 7 ) 

3=1 3=1 

where the subscript BC denotes the boundary conditions. In the transfer-matrix 
formalism, the Aqj's in the denominator of the generating function, eq. ([3.7)1 . are the 
eigenvalues of T. 

Strips of the triangular lattice are well-defined for widths L > 2. The partition 
function Z(G,q,v) has been calculated (for arbitrary q, v, and m) for the strip with 
L = 2p in 0] by two of us using a systematic iterative application of the deletion- 
contraction theorem. Z(G,q,v) was also studied for arbitrary q and v and zeros 
calculated for L = 3,4 and various lengths in jH]. Here after reviewing the Lp = 2p 
case in our formalism, we shall report explicit results for the partition function for 
Lp = 3f- For 4p < Lp < 6p, the expressions for T(Lp), w(Lp) and U;d(Lp) are so 
lengthy that we cannot include them here. They are available from the authors on 
request and in the MATHEMATICA file transf er_Tutte_tri .m which is available with 
the electronic version of this paper in the cond-mat archive at http://arXiv.org, 



3.1 L = 2 F 

Although this partition function was given earlier in jj], it is useful to review the 
calculation from the point of view of the (spin-representation) transfer matrix. The 
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number of elements in the basis is equal to C 2 = 2: P = {1, ^1,2}- In this basis, the 
transfer matrices and the other relevant quantities are given by 

/ q 2 + Aqv + 5v 2 + v 3 (1 + v) (q + 3v + v 2 ) \ 
' ~ V v 2 {q + 3v + v 2 ) v 2 {l + v)(2 + v) J 

w T = q{q + v, 1 + v) (3.8b) 

ul = (1,0) (3.8c) 

Because certain expressions recur in transfer matrices for wider strips, it is convenient 
to re-express ()3.8|) in terms of these expressions; we have 

T = ( Jk fo E k ) (Ma) 

w T = q(F x ,D x ) (3.9b) 

where 

D k = v + k (3.10a) 

F k = q + kv (3.10b) 

E k = q 2 + kv + v 2 (3.10c) 

T = q 2 + 4qv + 5v 2 + v 3 (3.10d) 

In terms of this transfer matrix and these vectors one calculates the partition function 
Z(G m ,q,v) for the strip with a given length m via eq. (J3.1|) . Equivalently, one can 
calculate the partition function using a generating function, and this was the way in 
which the results were presented in with 

2 

Af= JJ(l-A tAi *) (3.11) 

3=1 



and 

where 
and 



\ _ 1 

A 4 /2,(l,2) - - L 



T sl2 ±(q + 3v + v 2 )^/R^~ 2 \ (3.12) 
T S12 = v A + 4v 3 + 7v 2 + Aqv + q 2 (3.13) 



R S12 = q 2 + 2qv - 2qv 2 + 5v 2 + 2v 3 + v 4 (3.14) 
The product of these eigenvalues, which is the determinant of T, is 

det(T) = v 2 (l + v) 2 (v + q) 2 = v 2 D\Fl (3.15) 

The vanishing of this determinant at v — — 1 and v = —q occur because in each 
case one of the two eigenvalues is absent for, respectively, the chromatic and flow 
polynomials [37]. Analogous formulas can be given for det(T) for higher values of L; 
we omit these for brevity. 
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3.2 L = 3 F 



The number of elements in the basis is equal to C 3 = 5: P = {1, <5i j2 , #1,3, £2,3, ^1,2,3}- 
In this basis, the transfer matrices and the other relevant quantities are given by 



T 



w 



/ 



T\ 

v 2 F 3 E 3 

v 3 E 3 
v 2 E 3 E 4 
v 4 D 3 E 3 

<2 



V 



D\F 3 E 3 
2 D 1 D 2 F 2 
v 2 D 1 E 3 
v 2 D 1 D 2 E 3 
v 3 D 1 D 2 E 3 



U: 



id 



3 ^3 

q(F 2 ,D 1 F 1 ,E 2 ,D l F 1 ,Dl) 
1,0,0,0,0) 



T 2 
v 2 T 5 
v 2 T 5 
v 2 T 6 
v 3 T 7 



D\T 3 

• 2 D 1 T i 
,3, 



v 

v 3 D 1 D 2 



v 2 D 1 D 2 E. i 



D 2 T 4 

v 2 D 2 D 2 

v 2 D 2 D 2 

v 2 D 2 D 2 
3 r>2 n2 



l D t D 2 D 3 v A D[D 2 2 ) 



(3.16a) 

(3.16b) 
(3.16c) 



where the are shorthand notations used in this section, defined as 



Tl 
T 2 
T 3 
T 4 
T 5 
T 6 
T 7 



q 3 + 7q 2 v + 19qv 2 + 19v 3 + 2qv 3 + 5v 4 

q 2 + 7qv + 16v 2 + qv 2 + 9v 3 + 2v 4 

q 2 + 6qv + llv 2 + qv 2 + 4v 3 

q + 5v + 2v 2 

q + 6v + Av 2 + v 3 

2q + 13v + 13v 2 + 6v 3 + v 4 

q + I2v + 13f 2 + 6v 3 + v 4 



(3.17a) 
(3.17b) 
(3.17c) 
(3.17d) 
(3.17e) 
(3.17f) 
(3.17g) 



4 Potts Model Partition Functions for Triangular- 
lattice Strips with Cylindrical Boundary Condi- 
tions 

The Potts model partition function Z(G, q, v) for a triangular-lattice strip of width 
L vertices and length m vertices with cylindrical boundary conditions can be written 

as 

Z(L P x m F , q, v) = v T ■ H ■ T^ 1 • u id (4.1a) 
= w T -T m - 1 -u id (4.1b) 

where again w T = v T • H. Results on Z(G,q,v) for 2p < Lp < 4p with cylindrical 
boundary conditions were previously given in [5]. Here we shall present explicit results 
for this partition function for strips of the triangular lattice with cylindrical boundary 
conditions and widths 2p < Lp < 4p. Results for the transfer matrices in computer- 
readable files for widths up to Lp = 9p are also available upon request; they are too 
lengthy to present here. 

In the computation of the transfer matrix for a triangular-lattice strip with cylin- 
drical boundary conditions there is a technical complication in order to treat correctly 
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the last diagonal bond joining columns L and 1 jTHJ Section 3]. Instead of consider- 
ing a triangular-lattice strip of width L and cylindrical boundary conditions, we start 
with a strip of width L + l and free boundary conditions. The parameter v is the same 
for all edges, except for the vertical edges corresponding to column L + l, where it 
takes the value v — 0. After performing the computation, we identify columns 1 and 
L + l. This procedure implies that there are double vertical edges (with parameters 
v and 0) connecting sites on column 1; but their net contribution is v, as expected. 

We obtain in this way a transfer matrix of dimension Cl+i- This matrix can be 
simplified by noting that there are many zero eigenvalues. Let us denote by {v^ } 
(resp. {Vj } ) the basis elements corresponding to L + 1 being (resp. not being) a 

singleton. The number of elements of {v^} is Cl', hence the cardinality of {v^} 
is Cl+i ~ Cl- The zero eigenvalues are associated with certain eigenvectors of the 
form Vi — V2 where Vi G {vj } and V2 G {vj"' ) }. Let us make the following change 
of variable 

, fvf'-vf i = l,...,C L +i-C L 




where the first elements corresponds to eigenvectors with zero eigenvalues, and the 
last Cl elements are those of the basis {vj }. Then, the transfer matrix takes the 
simple form 

( J (n) \ 

T '=(°T(->) <«) 

where T^ (resp. T^>) is a matrix of dimension Cl (resp. Cl+i — Cl)- In this new 
basis the vectors w and take the form 

w' = (0,w (s) ) (4.4a) 
uL, = (0,ug } ) (4.4b) 

Thus, we can write the partition function in terms of T^, and alone 

Z(L P x m F , q, v) = w (s) T ■ [T^] m_1 ■ ug } (4.5) 

So far we have a transfer matrix with the same dimension as for free boundary 
conditions (namely, iVz itri ,FF,L = Cl)- We can reduce even more the dimension of the 
transfer matrix by noting that cylindrical boundary conditions introduce an extra 
symmetry (i.e. translations along the transverse direction). In particular, we can 
make a further change of basis in the subspace {v^} so that 



(s), = I v£'"> i = l,...,C L - N ZM>PF>L ( } 

^ vf' 4) i = C L -N ZM!PFi L + l,...,C L 



where the last Nz,tri,PF,L (resp. the first Cl — iVz,tri,PF,L) elements {v- s '*' ) } (resp. 
{v- s ' ra ' ) }) are translational-invariant (resp. non-translational-invariant) combinations 
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of the original vectors {v£ s '}. In this new basis, the transfer matrix takes a block 
diagonal form 

/ -r(s,n) n \ 

T< *"=( T<%) < 47 > 
and the vectors and take the form 

w (s) ' = (0,w (s,t) ) (4.8a) 

u!f = (0,u[f) (4.8b) 

The partition function can be computed using the transfer matrix T^'*^ 

Z(L P x m F , q, v) = T • [T^]" 1 " 1 • ujf } (4.9) 

The dimension of T^'*' will be denoted by Nz,tri,PF,L and is given by Theorem 12.21 fsee 
also Table HJ). 

For L > 6 there is a further simplification, which is also present in the chromatic- 
polynomial case [3UJ : in the translation-invariant subspace the transfer matrix does 
commute with the reflection operation. Thus, we can pass to a new basis consisting 
of connectivities that are either even or odd under reflection. In this new basis, the 
transfer matrix T^'*) has the block diagonal form 

*"* = ( T f T l> ) (4-10) 

and the vectors w^ s,t ^ and u^'^ take the form 

w^' = (0,w^ } ) (4.11a) 

vtfy = (0,nSa) (4.11b) 

Thus, we can compute the partition function using the transfer matrix T^'*)' 

Z(L P x m F , q, v) = wf > T • [T^]- 1 • u&J (4.12) 

In what follows, we shall drop the superindices ^ and the subindex + to simplify the 
notation. The dimension of the transfer matrix T+'^ = T corresponds to the number 
of partition classes of the numbers {1, 2, ... , L} that are invariant under translations 
and reflections. Thus, this number is precisely iVz iSqi PF,L and is given in terms of L 
by eq. ()2.52|) (see Tabled] for some numerical values). 

We have obtained the transfer matrices T(Lp) and the vectors w(Lp) and Ujd(Lp) 
using symbolic computation with mathematica for Lp < 5p. For 6p < Lp < 
9p, we have used two different programs (one written in C and the other one in 
perl), which were also used to double-check the results for L < 5p. For 5p < L < 
9p, the expressions for T(Lp) and the vectors w(Lp) and U;d(Lp) can be found in 
the MATHEMATICA file transf er_Tutte_tri . m which is available with the electronic 



version of this paper in the cond-mat archive at http://arXiv.org 
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4.1 L = 2 P 

The number of elements in the basis is two, with P = {1,5 12 }; the transfer 
matrices and the other relevant quantities are given by 

/ q 2 + 6qv + I2v 2 + qv 2 + 8v 3 + 2v 4 D 2 (q + 4t» + 2v 2 ) \ , , 
' ~ V v 2 (2q+12v+l3v 2 + 6v 3 + v i ) v 2 D 2 D 2 J {AA6a) 

w T = q(E 2 ,Df) (4.13b) 
"5 = (1,0) (4.13c) 

As before, we have used the shorthand notation introduced above ()3.10aj) / (|3.1Uc|) . 



4.2 L = 3 P 

The number of (translational-invariant) elements in the basis is three: P = 

{l,$i, 2 + $1,3 + 82,3, ^1,2,3}- The transfer matrices and the vectors w and u i( j are 
given by 

/ T n 3D^ l2 D 3 T 13 \ 

v 2 T 21 v 2 D 1 T 22 v 2 D\D\ (4.14a) 
\ v*D 3 T 31 3v 3 D 1 D 2 T 32 v 3 D 3 D 3 J 

q(V 1 ,3D 1 E 2 ,D 3 1 ) (4.14b) 

(1,0,0) (4.14c) 

where the factors D^, E& are given in (|3.1()ajl / (j3.10c|) . and 



Vi 


= q 2 + 3qv + 3v 2 + 2v 3 


(4.15a) 


T u 


= q 3 + 9q 2 v + 33qv 2 + Aqv 3 + 50v 3 + 2lv A + 3v 5 


(4.15b) 


T 12 


= q 2 + 8qv + 20v 2 + 2qv 2 + Uv 3 + 3v 4 


(4.15c) 


Tn 


= q + 6v + 3v 2 


(4.15d) 


T 2 i 


= q 2 + lOqv + 30v 2 + 2qv 2 + 22v 3 + 7v 4 + v 5 


(4.15e) 


T 22 


= 6q + 38v + 2qv + A2v 2 + 18v 3 + 3v A 


(4.15f) 


Tn 


= 3q + 18v + 15v 2 + 6v 3 + v 4 


(4.15g) 


T 32 


= q + 12v + 13v 2 + 6v 3 + v 4 


(4.15h) 



T 



w 



id 



4.3 L = 4 P 

The number of (translational-invariant) elements in the basis is six: P = {1, $\ t2 + 

$2,3 + $3,4 + $1,4, $1,3 + $2,4, $1,2,3 + $1,2,4 + $1,3,4 + $2,3,4, ^1,2,3,4, $1,4$2,3 + $1,2$3,4}- The 

transfer matrix is given by 

T = (4.16) 
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/ 
/ 




AD 1 1 12 




4:D 1 1 U 


-^1^15 


2L> 1 li e 




v 2 T 2l 


2v 2 D x T 22 


v 2 T 23 


v 2 D 2 T 24 


v 2 D 4 D 2 


2v 2 D 2 D 2 T 26 




v 3 T 31 


2v 2 D 1 T 32 


4^ 2 T| 4 


Av 2 D 2 D 2 T u 


v 2 D\D 2 


v 2 D 2 T 36 




v 4 D 3 T n 


v 3 D{T i2 


2v 3 T u T 54 


2v 3 D 2 D 2 T M 


v 3 D 4 D 3 


v 3 D 2 D 2 T m 




v 6 D 2 T 51 


Av^D x D 2 D 3 T^ 2 


2v% 2 4 


Av 4 D 2 D 2 T 54 


v 4 D 4 D 4 


2v 5 D 2 D 2 D 3 T 56 


V 


v 4 T 2 


Av 4 D 1 D 2 T m 


v 4 T 63 


2v 4 D 2 D 2 





2v 4 D\D\ 



he factors Dk{v) and Ek are defined in ()3.10a|) / ()3.10c|) ; the ly are given by 

T n = q 4 + I2q 3 v + 62q 2 v 2 + 16Aqv 3 + Aq 2 v 3 + 192w 4 + 29gv 4 

+72v 5 + 6v 6 (4.17a) 

T 12 = q 3 + llq 2 v + A7qv 2 + q 2 v 2 + 80v 3 + llqv 3 + AOv 4 + 5v 5 (4.17b) 

T 13 = q 2 + 8qv + 2Av 2 + qv 2 + 16v 3 + Av 4 (4.17c) 

T u = q 2 + Wqv + 28v 2 + 3qv 2 + 20v 3 + Av 4 (4.17d) 

T 15 = q + 8v + Av 2 (4.17e) 

T 16 = q 2 + lOqv + 32v 2 + 3qv 2 + 2Av 3 + 5v 4 (4.17f) 
T 21 = g 3 + 12q 2 v + 5Aqv 2 + 2q 2 v 2 + 96v 3 + 18qv 3 + 59v 4 + qv 4 

+12w 5 + w 6 (4.17g) 
T 22 = 3q 2 + 28qv + q 2 v + 80v 2 + 13qv 2 + 65v 3 + qv 3 + 18v 4 + 2v 5 (4.17h) 
T 23 = 3q 2 + 3Aqv + 96v 2 + 23qv 2 + 120v 3 + 8qv 3 + 67v 4 + qv 4 

+18v 5 + 2v 6 (4.17i) 

T 24 = lOq + 56v + 6qv + 63v 2 + qv 2 + 26v 3 + Av 4 (4.17j) 

T 26 = q + 8v + 5v 2 + v 3 (4.17k) 

T 31 = 2q 2 + I8qv + A8v 2 + Aqv 2 + 3lv 3 + 8t» 4 + v 5 (4.171) 

T 32 = q 2 + Uqv + 52v 2 + Aqv 2 + A3v 3 + Uv 4 + 2v 5 (4.17m) 

T 34 = g + 6w + Av 2 + v 3 (4.17n) 

T 36 = 2g + 24w + 27v 2 + 12v 3 + 2v 4 (4.17o) 

T 41 = g 2 + llqv + 36v 2 + 2qv 2 + 2Av 3 + 7v 4 + v 5 (4.17p) 
T 42 = 2g 2 + A3qv + q 2 v + 216v 2 + 28qv 2 + 270v 3 + Aqv 3 + 138v 4 

+36v 5 + Av 6 (4.17q) 

T 44 = 3g + 2Av + qv + 27v 2 + 12v 3 + 2v 4 (4.17r) 

T 46 = q + 2Av + 27v 2 + 12v 3 + 2v 4 (4.17s) 

T 51 = 4g + 2Av + 17v 2 + 6v 3 + v 4 (4.17t) 

T 52 = 2g + 18v + Ihv 2 + 6t» 3 + v 4 (4.17u) 

T 54 = q + 12v + 13v 2 + 6v 3 + v 4 (4.17v) 

T 56 = A + 3v + v 2 (4.17w) 

T 62 = g + 6w + 2w 2 (4.17x) 

T 63 = 3g + 24w + 26i; 2 + 12v 3 + 2v 4 (4.17y) 
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The vectors w and Ujd are given by 



w 



T 



q (V u mV 2 , 2E 2 2 , 4D 2 E 2 , Df, 2D 2 E 2 ) 
(1,0,0,0,0,0) 



(4.18a) 
(4.18b) 



u 



T 
•id 



where 



Vi 

v 2 



q 3 + Aq 2 v + 6qv 2 + 4f 3 + v 4 
q 2 + 3qv + 3v 2 + v 3 



(4.19a) 
(4.19b) 



5 Partition Function Zeros in the q Plane 

In this section we shall present results for zeros and continuous accumulation 
sets B q (in the g-plane) for the partition function of the Potts antiferromagnet on 
triangular-lattice strips of widths L < 5 with free or cylindrical boundary conditions. 

In Figure n w e show the partition-function zeros in the g-plane for strips of sizes 
Lp x (lOL)p with 2 < L < 5 and free boundary conditions. We also show the cor- 
responding limiting curves B q (L,v) for the infinite-length limit. Figure ^a) displays 
the zeros for v — — 1 (i.e., the chromatic zeros jlHUnUEOl), and Figures db)-(d) dis- 
play the corresponding zeros for the non-zero temperatures v = —0.75 (b), v = —0.5 
(c), and v = —0.25 (d). The case L = 2 was studied in |lj and the cases L = 3,4 
were studied in Figs. 3.15, 3.22] for arbitrary temperature, and our zeros are in 
agreement with these earlier works. 

The corresponding partition-function zeros and accumulation sets for triangular- 
lattice strips with cylindrical boundary conditions are shown in Figure El Again, 
Figure |2^a)-(d) displays the zeros for v = — 1 (a), v = —0.75 (b), v — —0.5 (c), and 
v = —0.25 (d). The cases L = 3,4 were studied previously in (Figs. 3.19, 3.25 of) j3] 
for arbitrary temperature, and our zeros are in agreement with this work. 

The case v — — 1, which is the zero-temperature Potts antiferromagnet (chromatic 
polynomial) has been previously studied in [1HJ HOI EQ 113 HiO] for the free longitudi- 
nal boundary conditions and in ESI EH EH 120] for periodic longitudinal boundary 
conditions. For the case of free longitudinal boundary conditions [SO] contains results 
for L < 9f and L < 12p. Our Figures Ufa) andEfa) include calculations up to L = 5 
for comparison with other values of v. Although some curves, such as those for cylin- 
drical boundary conditions, may enclose regions, the curves do not enclose regions 
containing the origin. One observes that for either type of transverse boundary con- 
dition, as the width L increases, the left-hand arc endpoints move slowly toward the 
origin. When this was observed in earlier work for several different lattice strips, [48j, 
it motivated the suggestion that in the limit L — > oo for strips with free longitudinal 
boundary conditions, the limiting B q would separate the q plane into regions includ- 
ing a curve passing through q = jlHJ [T7j. The specific calculation of B q {v = — 1) 
in the limit L —>■ oo reported by Baxter j^B] has this feature. (For critical comments 
on certain features of Baxter's results, see the detailed discussion in [30J.) The prop- 
erty that B q separates the q plane into regions with one of the curves on B q passing 
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through the origin is also observed for lattice strips with finite width L if one imposes 
periodic longitudinal boundary conditions |H3 EH1 E21 EH1 EH EE Ell EH EH] • 

In making inferences about possible L — > oo characteristics of the continuous 
accumulation set of zeros B(G S , L x oo) in the q or v plane for infinite- length, width- L 
lattice strip graphs of type G s , one should recall that, in general, lim^oo B(G S , L x 
oo) is different from the continuous accumulation set of the zeros of the partition 
function for the usual 2D thermodynamic limit defined by starting with an L y x 
L x section of a regular lattice and letting L x and L y both approach infinity with 
L y /L x a finite nonzero constant. This type of noncommutativity was encountered 
in previous studies of B for the Potts model free energy on infinite-length, finite- 
width strips with periodic longitudinal boundary conditions [THJ EJ I2H1 122] ; for these 
strips, B is noncompact in the v plane, reflecting the fact that the Potts model has a 
ferromagnetic critical point only at T = (i.e. K = oo, hence v = oo) for any width 
L, no matter how great, whereas for the 2D lattice defined in the thermodynamic 
limit, it has a ferromagnetic critical point at a finite temperature, so B is compact 
in the v plane. Noncommutativity of this type was also found in studies of B q . 
For example, in calculations of B q for infinite-length strips of the triangular lattice 
with cyclic boundary conditions, it was found that this locus always passes through 
q = 2 [521 E31 EH El], whereas, in contrast, the locus found in [46J for the infinite- 
width limit of strips with cylindrical boundary conditions does not pass through 
q = 2. Similarly, in calculations of B q for infinite-length strips of the square lattice 
with cyclic boundary conditions, it was found that this locus always passes through 
q = 2, whereas in calculations of B q for infinite-length strips of the square lattice 
with cylindrical boundary conditions j5UJ HH1 EH1 E2j , it was found that B does not 
pass through q = 2, strongly suggesting that this difference will persist in the limit 

Ly — > OO. 

At nonzero temperature for the antiferromagnet, as represented in our figures by 
the range —0.75 < v < —0.25, the partition-function zeros in the g-plane have a 
different shape: as L increases, the limiting curves B q (L,v) tend to a bean-shaped 
curve or set of arcs, open on the left, without substantial protruding branches, in 
contrast to many of the v = —1 curves. For a given value of v in the range considered, 
as the width L increases, the curve envelope moves outward somewhat and the arc 
endpoints on the left move slowly toward q = 0. This behavior is consistent with 
the hypotheses that for a given v , as L — > oo, (i) B q would approach a limiting 
locus as L — > oo and (ii) this locus would separate the q plane into different regions, 
with a curve passing through q = as well as a maximal real value, q c {v). This is 
qualitatively the same type of behavior that was found earlier for the square-lattice 
strips [T^l 1231 13] - I n particular, our results are consistent with the inference that as 
L — > oo, B q for v = — 1 would pass through q c (v = —1) = 4, corresponding to the 
property that the q = 4 Potts antiferromagnet has a zero-temperature critical point 
on the (infinite) triangular lattice [60J. 

For a given L, as v increases from —1 to 0, i.e., as the inverse temperature f3 
decreases from infinity to for the antiferromagnet, the zeros and the limiting curve 
contract to a point at q — 0. This is an elementary consequence of the fact that 
these lattice strips have fixed maximal vertex degree and as the parameter K = f3J 
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approaches zero, the spin-spin interaction term in TC, eq. (jl.2|) . vanishes, so that the 
sum over states just counts all q possible spin states independently at each vertex, and 
Z(G, q, v) approaches the value Z(G, q, 0) = q n . An upper bound on the magnitudes 
of the zeros is given by the following theorem by Sokal: 

Theorem 5.1 (Sokal L 611) Let G = {V, E} be a loopless finite undirected graph of 
maximum degree < r, equipped with complex edge weights {v e } eG E satisfying \ l + v e \ < 
1 for all e. Let |t>| max = max eg £ |t> e | . Then all the zeros of Z(G,q, {v e }) lie in the 
disc \q\ < C(r> max with C(r) < 7.963907r. 

This is a loose bound; for all of the strips with cylindrical boundary conditions 
and for the free strips with widths L > 3, the maximal degree is r = 6, so that the 
above theorem implies that |g| < 47.8|t> | max . Thus, for example, for v = —1, this 
reads \q\ < 47.8, whereas in fact \q\ < 4 for free boundary conditions and |g| < 4.5 
for cylindrical boundary conditions. A general feature is that the limiting curves and 
associated zeros tend to be located mostly in the Re(g) > half plane. 

As noted, it is evident in Figures ^ and 121 that as L increases, the accumulation set 
B q (L, v) moves outwards. As expected, the convergence to the limit L —>■ oo seems to 
be faster with cylindrical boundary conditions, as there are no surface effects when 
the length is made infinite. 

One can also plot B q for the ferromagnetic region < v < oo (e.g., Figs. 2 and 13 
of jlj). Although we have not included these plots here, we note that an elementary 
Peierls argument shows that the g-state Potts ferromagnet on infinite-length, finite- 
width strips has no finite-temperature phase transition and associated magnetic long 
range order. Hence, for this model B q does not cross the positive real q axis for 
< v < oo. 

6 Partition Function Zeros in the v Plane 

6.1 General 

In this section we shall present results for zeros and continuous accumulation 
sets B v (in the w-plane) for the partition function of the Potts antiferromagnet on 
triangular-lattice strips of widths L < 5 and free or cylindrical boundary conditions. 
Our results hold for arbitrarily great length L and for any real or complex value of 
q; they thus complement calculations of the Potts model partition function for fixed 
positive integer values of q on sections of the triangular lattice jH21 ESI EH ESI- We 
shall focus here on integer values of q, since these are the most relevant from a physical 
point of view. We recall the possible noncommutativity in the definition of the free 
energy for certain integer values of q (see eqs. (2.10), (2.11) of [TH] or (1.17) of jl]): 
lim^oo \im q __> qs Z(G, q, v) x l n ^ Yim. q _ >qa Yim n ^ 00 Z{G,q,v) l / n As discussed in |T8] . be- 
cause of this noncommutativity, the formal definition (jl.4j) is, in general, insufficient to 
define the free energy / at these special points q s ; it is necessary to specify the order of 
the limits that one uses in the above equation. We denote the two definitions using dif- 
ferent orders of limits as f qn and f nq : f nq ({G}, q, v) = lim^oo lim^^ n^ 1 In Z(G, q, v) 
and f qn ({G},q,v) = lim g _> 9s lining n _1 In Z(G, q, v). 
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As a consequence of this noncommutativity, it follows that for the special set 
of points q = q s one must distinguish between (i) (B v ({G},q s )) nq , the continuous 
accumulation set of the zeros of Z(G,q,v) obtained by first setting q = q s and then 
taking n — > oo, and (ii) (B V ({G}, q s )) qn , the continuous accumulation set of the zeros 
of Z(G,q,v) obtained by first taking n — > oo, and then taking q — ► q s . For these 
special points (cf. eq. (2.12) of [T5]). 

(B V ({G}, q s ))nq ^ (B V ({G}, q s )) qn (6.1) 

Here this noncommutativity will be relevant for q = and q = 1. 

In Figures OHZI we show the partition-functions zeros in the v -plane (for a fixed 
value of q) for strips of sizes Lp x (10L)p with 2 < L < 5 and free boundary conditions. 
We also show the corresponding limiting curves B V (L, q) for the limit of infinite strip 
length. For simplicity, we have displayed each value of L on a different plot: L — 2 (a), 
L — 3 (b), L — 4 (c), and L = 5 (d). The corresponding partition-function zeros and 
accumulation sets for triangular-lattice strips with cylindrical boundary conditions are 
shown in Figures IH1fT2l with the same notation as for the former figures. Complex- 
temperature phase diagrams and associated partition function zeros were given in ^] 
for L = 2 for free and periodic longitudinal boundary conditions and free transverse 
boundary conditions. Results for L = 4p and L = 4p at zero and finite temperatures 
were given before in [3] . Our present calculations are in agreement with, and extend, 
this previous work. 

On the infinite triangular lattice (defined via the 2D thermodynamic limit as given 
above), the phase transition point separating the paramagnetic (PM) and ferromag- 
netic (FM) phases is determined as the (unique) real positive solution of the equation 

m 

v 3 + 3v 2 -q = (6.2) 

In previous studies such as [T%1 H] , it has been found that although infinite-length, 
finite-width strips are quasi-one- dimensional systems and hence the Potts model 
has no physical finite-temperature transition on such systems, some aspects of the 
complex-temperature phase diagram have close connections with those on the (infi- 
nite) triangular lattice. We shall discuss some of these connections below. 

6.2 q = 

From the cluster representation of Z(G, q, v), eq. (jl.5|) . it follows that this partition 
function has an overall factor of q k<yG \ where k(G) denotes the number of components 
of G, i.e., an overall factor of q for a connected graph. Hence, Z(G,q — 0, v) = 0. 
In the transfer matrix formalism, this is evident from the overall factor of q coming 
from the vector w. However, if we first take the limit n —>■ oo to define B for q ^ 
and then let q — > or, equivalently, extract the factor q from the left vector w, 
we obtain a nontrivial locus, namely (B V ({G}, 0) qn . This is a consequence of the 
noncommutativity ()6.1|) for q = 0. 

With the second order of limits or the equivalent removal of the factor of q in 
Z, we obtain the locus B v (q = 0) shown in Figures El (free boundary conditions) 
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and |H1 (cylindrical boundary conditions). The accumulation set B v (L,q = 0) seems 
to converge to a roughly circular curve. We see in Figures H3 and |H1 that the limiting 
curves cross the real v-axis at v ~ —3. We note the interesting feature that this 
is a root of eq. (|6.2j) for q = 0. For the case of cylindrical boundary conditions, 
B v (q = 0) includes a small line segment on the real axis near v = —3, with a length 
that decreases as L increases. As L increases, the arc endpoints on the upper and 
lower right move toward the real axis. It is possible that these could pinch this axis 
at v = as L — > oo, corresponding to the other root of (|6.2j) for q = 0. 

6.3 q = l 

For q = 1, the spin-spin interaction in (jl.2j) always has the Kronecker delta func- 
tion equal to unity, and hence the Potts model partition function is trivially given 
by 

Z(G,q = l,v) = e K \ R \ = {l + v)W (6.3) 

where \E\ is the number of edges in the graph G. This has a single zero at v — — 1. 
But again, one encounters the noncommutativity (J6.1|) for q — 1. It is interesting to 
analyze this in terms of the transfer matrix formalism. At this value of q, both the 
transfer matrix and the left vector w are non-trivial. There is thus a cancellation of 
terms that yields the result (j6.3|) . The case L = 2p is the simplest one to analyze: 
the eigenvalues and coefficients for q = 1 are given by 

Ai(M) = 2v 2 ; Cl (l,0) = (6.4a) 

A 2 (M) = (l + v) 6 ; c 2 (l,0) = (l + i;) 2 (6.4b) 

Thus, only the second eigenvalue contributes to the partition function, and it gives 
the expected result Z(2 P x m F , q = 1, v) = (1 + t>) 6m_4 . For L = 2 F we obtain: 

Ai(l,v) = v 2 ; ci(l,0) = (6.5a) 

\ 2 {l,v) = (l + v) 4 ; c 2 (l,0) = (1 + v) (6.5b) 

giving rise to Z(2p x trip, q = 1, v) = (1 + v ) 4m_3 . The case L = 3p is similar: there is 
a single eigenvalue Ai(l,t>) = (1 + v) 9 with a non-zero coefficient c\(l,v) = (1 + f) 3 , 
and the other two eigenvalues [which are the roots of x 2 — xv 2 (3v 3 + 13v 2 + 24v + 
3) + v 5 (1 + f )(6 + 5v + 6v 2 )} have identically zero coefficients 02,3(1, v) = 0. Thus, the 
partition function takes the form Z(3p x mp, 9 = 1, = (1 + v) 9m ~ 6 = (1 + v )' s '. 

In general, we conclude that at q — 1 only the eigenvalue A = (1 + v ) 3Lp (resp. 
A = (1 + f ) 3Lp ~ 2 ) contributes to the partition function for cylindrical (resp. free) 
boundary conditions, and its coefficient is c = (1 + v) Lp (resp. c = (1 + t>) iF_1 ). 
The other eigenvalues do not contribute as they have zero coefficients. This is the 
analogue of what was found for the strips with cyclic boundary conditions, where the 
various A's fall into sets \G,dj such that all of the A's with a fixed d have a unique 
coefficient which is a polynomial of degree d in q given by [T%1 |20] 
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where U u (z) is the Chebyshev polynomial of the second kind. These coefficients 
vanish at certain values of q, which means that if one evaluates the partition function 
first at these values and then takes the limit n — > oo, the corresponding A's will 
not contribute to Z, while if one takes n — > oo first, calculates the free energy and 
the locus B qn , and then sets q equal to one of these values, the A's will, in general 
contribute. In particular, we recall (eq. (2.18) of [201) that ^ Q = 1> then c^ d > vanishes 
if d = 1 mod 3. Thus, we see similar manifestations of the noncommutativity (j6.1|) 
for strips with free and periodic boundary conditions. 

In our present case, in order to obtain B qn , we have computed the Tutte-polynomial 
zeros and the corresponding limiting curves for q = 0.999 (see Figures 0] and EJ) • The 
accumulation sets B v (L,q = 0.999) for L = 2 F to L = 5p consists of arcs that come 
close to forming an almost closed bean-shaped curve, with an involution on the right 
that deepens as L increases. Solving the q = 1 special case of eq. (|6.2j) yields the 
roots v = -2.879385.., v = -0.652703..., and v = 0.5320888... The locus B v crosses 
the real v axis at two points, and our results are consistent with the inference that 
as L — > oo, these two crossing points are the first two roots listed above. For free 
boundary conditions with width L = 3,4,5, B qn exhibits a small involution on the 
left. A noteworthy feature of this locus is that it is relatively smooth, without the 
prongs that tend to occur for the other values of q discussed here. 

6.4 q = 2 

The zeros and accumulation sets for q = 2 are displayed in Figures ISl and ITU1 for 
free and cylindrical boundary conditions. Fig. EJa) contains the same information 
as Fig. 4 of j3] (which is plotted in a different temperature variable, a -1 ). The 
finite-size effects for the accumulation sets B v (L,q = 2) are noticeably larger for 
free, in comparison with cylindrical, boundary conditions, as expected. In the latter 
case, the curves B v (L,q = 2) for L = 2,3,4,5 fall very approximately one on top 
of the preceding one. For cylindrical boundary conditions, we see that the curve 
B v (L,q = 2) is symmetric under the replacement a — > —a. The reason for this is 
that in this case all of the vertices except for the end-vertices, which constitute a 
vanishingly small fraction in the limit of infinite length, are equivalent (i.e., the graph 
is r-regular) and have even degree r. In general, this property applies for the complex- 
temperature phase diagram of the q = 2 (Ising) special case of the Potts model for 
an infinite lattice where the coordination number is even jHHl ESI IZOj • Our strips with 
free transverse boundary conditions are not r-regular graphs because the vertices on 
the upper and lower sides have a different degree than those in the interior. Because 
of this, the B v in this case does not have the a — > —a symmetry. From previous 
work fSJ H] one knows that the loci B v are different for strips with free or periodic 
transverse boundary conditions and free longitudinal boundary conditions, on the one 
hand, and free or periodic transverse boundary conditions and periodic (or twisted 
periodic) longitudinal boundary conditions. One anticipates, however, that in the 
limit of infinite width, the subset of the complex-temperature phase diagram that 
is relevant to real physical thermodynamics will be independent of the boundary 
conditions used to obtain the 2D thermodynamic limit. 
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In the 2D thermodynamic limit, one knows the complex-temperature phase di- 
agram exactly for the q = 2 (Ising) case. (This isomorphism involves the redefini- 
tion of the spin-spin exchange constant Jp tts = 2Ji sing and hence Kp oUa = 2Kj sing , 
where i^p tts is denoted simply K here.) The simplest way to portray the complex- 
temperature phase diagram is in the a 2 or u = a -2 plane since this automatically 
incorporates the a — > —a symmetry noted above. In the u plane, the complex- 
temperature phase diagram, with boundaries given by B u , is (see Fig. 1(a) of [7U] 
which is equivalent, by duality to the complex-temperature phase diagram for the 
honeycomb lattice given as Fig. 2 in j72j ) 

B u : {| M + i| = ^} U {-oo<u<-^} (6.7) 

i.e., the union of a circle centered at u = — 1/3 with radius 2/3 and the semi-infinite 
line segment extending leftward from u = —1/3 along the real u axis. In the a -1 
plane (Fig. 1(b) of (70], related by duality to Fig. 3 of [72]), B is the union of a 
vertically elongated oval crossing the real axis at ±1/V3, the imaginary a -1 axis at 
±i, and two semi-infinite line segments extending from ij\f% to zoo and from 
to —ioo along the imaginary axis. Equivalently, in the a plane, B is the union of 
a horizontally elongated oval crossing the real a axis at ±a/3 and the imaginary a 
axis at ±i, and a line segment along the imaginary axis extending between \/3i and 
— \/3i. The locus B v in the v plane is obtained from this by translation by one unit, 
since v = a — 1. This locus separates the complex v plane into three phases: (i) 
the paramagnetic phase, including the infinite-temperature point v = 0, where the 
S q symmetry is realized explicitly (S q being the symmetric group on q numbers, the 
symmetry group of the Hamiltonian) , (ii) the ferromagnetic phase, including the real 
interval v c (q = 2) < v < oo where the S q symmetry is spontaneously broken by the 
existence of a nonzero magnetization, and (iii) an unphysical phase (denoted "O" for 
"other" in [70]) including the point v = — 2. Here 

a c (q = 2) = v c (q = 2) + 1 = V% (6.8) 

is the physical critical point separating the PM and FM phases (for a review of the 
Ising model on the triangular lattice, see, e.g., [73]) • These physical PM and FM 
phases have complex-temperature extensions off the real v axis. The PM and O 
phases are separated by the subset of the vertical line segment extending between 
a = i and a = —i; this line segment terminates at the points a = ±\^3i. Because of 
the maximal frustration, there is no antiferromagnetic phase at finite temperature. 
The presence of a zero-temperature critical point in the 2D Ising antiferromagnet [73] 
is manifested by the fact that B v passes through v = —1, i.e., a = (as part of the 
above-mentioned vertical line segment). The complex-temperature phase boundary 
B v crosses the real v axis at v — a/3 — 1, separating the FM and PM phases, at 
v — — 1, separating the PM and O phases, and at v — — 1 — \/3, separating the O 
and (complex-temperature analytic continuation of the) FM phases. In [3] the B for 
an infinite-length free or cyclic strip with width L = 2 were compared with this 2D 
phase diagram. These three points, v = — 1, — 1 ± a/3, are the three roots of the q = 2 
special case of eq. (|6.2jl . 
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Using our exact results, we can compare our loci B v for a wide variety of widths 
and either free or periodic transverse boundary conditions with the known complex- 
temperature phase diagram for the Ising model on the infinite 2D triangular lattice. 
This comparison is simplest for the case of cylindrical boundary conditions, so we 
concentrate on these results. For the finite values of L that we have considered, B v 
has the form of two complex-conjugate arcs that cross two complex-conjugate line 
segments on the imaginary axis at v — — 1 ± i. One sees that as L increases, the 
endpoints of the arcs move down toward the real axis, as do the endpoints of the 
line segments. As L — ► oo, we expect that these arc endpoints will close, forming 
the above-mentioned horizontally elongated oval and vertical line segment extending 
from v = — 1 + v3i to v = — 1 — y/3i that constitute the complex-temperature phase 
boundaries B for the Ising model on the infinite triangular lattice. 

6.5 q = 3 

In contrast to the q = 2 case, the free energy of the g-state Potts model has not 
been calculated exactly for q > 3 on any 2D (or higher-dimensional) lattice and hence 
the complex-temperature phase diagrams are not known exactly. The q = 3 special 
case of eq. (|6.2J) has the root 

vpM-FM, g =3 = v c (q = 3) = -1 + cos(2tt/9) + v / 3sin(2 7 r/9) = 0.879385... (6.9) 

corresponding to the physical PM-FM phase transition point, and two other roots at 
the complex-temperature values 

v = -l + cos(2tt/9) - V3sin(27r/9) = -1.347296... (6.10) 

and 

v = -l-2 cos(2tt/9) = -2.532089... (6.11) 

Discussions of the complex-temperature solutions of eq. ()6.2|) and their connections 
with the complex-temperature phase diagram have been given in IMl . A 

number of studies involving exact calculation of the partition function for various q 
values on large sections of the triangular lattice have been performed jfi^ l6~4l 163] . 
(There have also been many studies calculating zeros for the Potts model with q > 3 
on the square lattice; see [3] for references to these works.) 

The zeros and accumulation sets for q = 3 are displayed in FigureslfilandlTTlfor free 
and cylindrical boundary conditions. We expect that the pair of complex-conjugate 
endpoints in this regime will eventually converge to the ferromagnetic critical point 
v c(q — 3) as L — *■ oo. However, obviously, an infinite- length strip of finite width L is a 
quasi-one-dimensional system, so the Potts model has no physical finite-temperature 
phase transition on such a strip for any finite L. 

In the antiferromagnetic regime — 1 < v < 0, we observe noticeable finite-size 
effects even with cylindrical boundary conditions. In this regime, we also observe a 
complex-conjugate pair of endpoints with small value of lm(v) that, as L — > oo, are 
expected to approach the real v axis at the transition point separating the param- 
agnetic and antiferromagnetic (AFM) phases of the q = 3 Potts antiferromagnet on 
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the infinite triangular lattice. Monte Carlo and series analyses |75| [T^J [77J |7Hj have 
yielded the conclusion that the PM-AFM transition in the q = 3 Potts antiferromag- 
net on the triangular lattice is weakly first-order. A high-accuracy determination of 
the location of the PM-AFM transition temperature T was obtained in [7Hj by means 
of Monte Carlo simulations: T = 0.62731 ± 0.00006, or equivalently 

vp M -AFM,q=3 = -0.79691 ± 0.00003 (6.12) 

We shall improve this estimate below. 

Finally, in the complex-temperature interval v < — 1, the finite-size and bound- 
ary condition effects are evidently very strong. Because of this, in previous work, a 
combination of partition-function zeros and analyses of low-temperature series expan- 
sions was used jHS]; these enable one at least to locate some points on the complex- 
temperature phase boundary. As regards the infinite 2D triangular lattice, because 
of a duality relation, the complete physical temperature interval < T < oo, i.e., 
< a < 1 of the g-state Potts antiferromagnet on the honeycomb lattice is mapped 
to the complex-temperature interval — oo < v < — q) on the triangular lattice (and 
vice versa) |S5| . As was noted in f63., it follows that because the q = 3 Potts an- 
tiferromagnet on the honeycomb lattice is disordered for all temperatures, including 
T = 0, the free energy for this model on the triangular lattice is analytic in the inter- 
val — oo < v < —3, and hence no part of the complex-temperature phase boundary B v 
can cross the negative real axis in this interval. In particular, one anticipates that as 
L — > oo for the infinite-length, width- L strips, the left-most arcs on B v will not close 
and pinch the negative real axis in this interval. Our calculations of B v are consistent 
with the inference that as L — > oo, this locus crosses the real axis at the points (|6.10|) 
and (J6.ll)) . although there are significant differences between the loci obtained with 
free and cylindrical boundary conditions. 

We also observe certain line segments on the real v axis in the complex-temperature 
region. We note that massless phases with algebraic decays of correlation functions 
have been suggested for the Potts model on the (infinite) square lattice at real val- 
ues of v and q in the intervals —2 — y/A — q < v < — 2 + \/A — q with q G (0, 4) and 
q 7^ B r = A cos 2 (7r/r), and it was conjectured that these might also occur for other 2D 
lattices jM] • However, the correspondence of these suggestions with our results is not 
clear; for example, the above interval suggested in [HE] shrinks to zero as q — > 4, but 
we observe clear line segments on the real v axis for q = A for both free and cylindrical 
boundary conditions (see Figs. |7|and[T2J). A possible physical subset of the above 
range of v given in [66] would be the antiferromagnet interval — 1 < v < 0. However, 
the condition that q ^ B r excludes all of the integral values of q in the indicated range 
(recall that B2 — 0, B3 — 1, B4 — 2, B$ — 3, and B\ = lim n ^ 00 .B ri = 4). The claim 
in jnni is thus complicated by the fact that although it is possible formally to define 
the Potts model partition function Z(G, q, v) using (|1.5j) for real positive non-integral 
q for the antiferromagnetic case, — 1 < v < 0, here the model does not satisfy the 
usual statistical mechanical requirement that the partition function is positive, and 
hence does not, in general, admit a Gibbs measure (TTJUHj. This leads to pathologies 
that preclude a physical interpretation, such as negative partition function, negative 
specific heat, and non-existence of a \ V\ — > 00 limit for thermodynamic functions that 
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is independent of boundary conditions [7TJ 0j. As regards the connection with 
the locus B, a signal of a massless phase would be a line segment on B on the real v 
axis for fixed q or the real q axis for fixed v. For the zero-temperature Potts antifer- 
romagnet, i.e., chromatic polynomial, v = — 1, these phases would thus occur in the 
intervals between the Beraha numbers, < 5 < 1, l<g<2, 2 < q < (l/2)(3 + a/5), 
and so forth. However, it has been proved that there are no real zeros of a chromatic 
polynomial in the intervals —00 < q < 0, < q < 1, and 1 < q < 32/27 [7^1 IHUj . 
Since B forms as an accumulation set of zeros, this makes it difficult to see how there 
could be a line segments in these intervals, in particular, the intervals < q < 1 
and 1 < q < 32/27. Again, it is not clear how to reconcile the results of these theo- 
rems with a suggestion that there would be massless phases with associated real line 
segments on B in these intervals. 

6.6 q = A 

For q = 4, eq. (j6.2j) has the physical root 

vpM-FM, q =4 = v c (q = 4) = 1 (6.13) 

corresponding to the PM-FM phase transition point and a double root at the complex- 
temperature point 

v = -2 (6.14) 

The zeros and accumulation sets for q = 4 are displayed in Figures [7| and One 
observes the approach of the right-most complex-conjugate arcs to the real 
increases, i.e. the approach to the PM-FM critical point in this case. For a given L, 
the approach to the exactly known value vpM-FM, q =i = 1 in eq. (|6.13|) is closer for 
cylindrical versus free boundary conditions, as is anticipated since the former minimize 
boundary effects. The q = 4 Potts antiferromagnet on the triangular lattice has a zero- 
temperature critical point, so that v = — 1 is on B v [SU] (this is not a root of eq. ()6.2I) ). 
For L > 2 for free boundary conditions and for L > 2 for cylindrical boundary 
conditions, we see how a pair of complex-conjugate arc endpoints approaches the 
real axis as L increases, consistent with the inference that these would pinch at 
v = —1 in the L — > oo limit. For both free and cylindrical boundary conditions and 
various values of L, one sees that B v contains an intersection point at the complex- 
temperature value v = —2, in agreement with the expectation from eq. (|6.14j) . The 
fact that the q = 4 Potts antiferromagnet is disordered on the honeycomb lattice for 
all temperatures T including T = implies that B v does not cross the negative real 
axis in the interval — oo < v < —4 |63j . In particular, this implies that as L —>■ oo, the 
leftmost arc endpoints on B v in the figures do not move down to pinch the negative 
real axis in this interval — oo < v < —4, provided that the limit L —>■ oo of these 
infinite-length, width-L strips commutes with the 2D thermodynamic limit for the 
triangular lattice as regards this aspect of the complex-temperature phase diagram. 

We do not show plots for q > 5, but recall that the Potts antiferromagnet is 
expected to be disordered (with exponential decay of spin-spin correlation functions) 
even at T = on the triangular lattice. This can be proved rigorously for q > 11 as 
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a slight improvement of the result that g-state Potts antiferromagnet is disordered at 
all temperatures on a lattice with coordination number r if q > 2r [71]. The property 
that the Potts antiferromagnet is disordered at all T on the triangular lattice for 
q > 5 is reflected in the property that B v does not pass through v = —1. 

7 Internal Energy and Specific Heat 

The partition function (jl.ljl can be used to derive the free-energy density f(G, q, v) 

f(G,q,v) = -L\ ogZ (G,q,v) (7.1) 

for finite \V\, with the V — > oo limit having been defined in eq. (jl.4|) above. The 
internal energy E and the specific heat C are derived in the usual way from the free 
energy as 

E(G,q,v) = -J^ = -J(v + l)^ (7.2) 

and 



dE , ^ \df , , , lN <9 2 / 



-J- + [v + l) 



k B K 2 C H (7.3) 



(9f cfa; 2 

Henceforth, for convenience, we shall use a definition of E without the factor —J in 
(J7.2|) . and we shall use the dimensionless function Ch in discussions of the specific 
heat. Let us suppose that G is a triangular-lattice strip graph of size L x m. In 
the limit m — > oo, since only the dominant eigenvalue Xd{q,v) of the transfer matrix 
contributes to the free energy, one has 

= ^logA, (7.4a) 

£(<Z,<y;L) = ^Ai (7.4b) 

C H (q,v;L) = ^[A 2 + A 1 -Af\ (7.4c) 



where the Aj(g, f ; L) are defined by 



Ai{q, u; L) = — (7-5) 

The infinite-length limits of the triangular-strips considered here are quasi-one- 
dimensional systems with analytic free energies at all temperatures. Hence the dom- 
inant eigenvalue Ad is the same on the whole semi-axis Im(t>) = 0, Re(t>) > — 1. 
Furthermore, as discussed in the free energy and its derivatives with respect to 
the temperature are independent of the longitudinal boundary conditions in the limit 
m — > oo (although they depend on the transverse boundary conditions). 

In Figure HH1 we have plotted the internal energy E ()7.2|) , the specific heat Cjj 
(|7.3|) . and the Binder cumulant £7 4 (|7.11|) (see below) for q = 3 on a triangular-lattice 
strip of width up to L = 6 and infinite length with cylindrical boundary conditions. 
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The behavior of the energy for the triangular-lattice strips with cylindrical bound- 
ary conditions is interesting: the curves cross each others close to the critical value v c 
in the ferromagnetic regime. In particular, for q = 2 we find that all curves cross at 
v c (2) = a/3 — 1 ~ 0.7320508.... For q = 3, 4 we find that the crossings are close to the 
respective PM-FM critical points t> c (3) given exactly in eq. (j6.9J) (= 0.8793852...) and 
f c (4) = 1, but they do not coincide precisely with these critical points. The reason for 
the above behavior is the following: the triangular-lattice Potts model on a triangular 
lattice at a given value of the temperature Boltzmann variable v is related by duality 
to the hexagonal-lattice Potts model at a different temperature variable 

v'=l (7.6) 
v 

This relation is exact when the original triangular-lattice is defined on an infinitely 
long cylinder of width L. In the Ising case q = 2, there is an additional transformation 
(namely, the star-triangle transformation) that maps back the hexagonal-lattice Potts 
model onto a triangular-lattice Potts model at temperature variable v"(v). This 
transformation allows us to compute the critical temperature v c (i.e., v c is the unique 
fixed point of the equation v"(v) = v) and the critical energy E c = E(q, v c ; L). When 
we perform the computation of E c in the limit m — ► oo, all finite-size corrections 
disappear, so that [HI] 

E c (q = 2-L) = E{q = 2, v c ; L) = E c {q = 2; oo) = - (7.7) 

However, for q ^ 2 there is no star-triangle transformation, and this implies the 
existence of corrections to scaling: 

oo . 

E c (q; L) = E(q, v c] L) = E c (q; 00) + ^-^, q ^ 2 (7.8) 

k=l 

where the parameters {uk(q)} are correction-to-scaling exponents that depend in 
general on the value of q. 

One can obtain a pseudo-critical temperature ve = VE(q,L,L') by solving the 
equation 

E(q,v E ,L) = E(q,v E ,L') (7.9) 

When L, L' — >• 00, we expect that this quantity will converge to the true critical value 
v c (q). This method has been employed in the the literature to locate critical points 
for several statistical- mechanical systems [S21 El and references therein] . 

Another pseudo-critical temperature can be obtained by looking at the point Vc = 
vc{q, L) where the specific heat C(q, v; L) attains a maximum value. This value differs 
from the bulk critical value v c (q) by finite-size-corrections of order ~ L~ x l v [HI] . 

We can also consider higher derivatives of the free energy with respect to K. In 
particular, the quantity M4 is the fourth derivative of the free energy with respect to 
K and, in the limit m — > 00 it can be written as 

u 4 (q,v;L)^-^ (7.10) 
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Instead of u^, it is more useful to deal with the phenomenological quantity U4 (also 
called the Binder cumulant [85]) defined as 

U^v;L) = ^ 2 (7.11) 

(Since this involves the energy instead of the magnetization, one could also consider 
a third cumulant.) A plot of this quantity for q = 3 is given in Figure ^fc). A third 
pseudo-critical temperature vu = vu{q,L) can be defined as the value at which the 
Binder cumulant (J7.11j) attains a minimum value. Again, this estimate is expect to 
differ from the bulk critical temperature v c (q) by terms of order L^ 1 ^ [S3]. 

The computation of the transfer matrices for triangular-lattice strips of width L 
allows us to compute these three pseudo-critical temperatures (namely, V£(q, L, L'), 
vc{q, L), and vu(q, L)). We have computed these estimates for several values of 1/2 < 
q < 4 in the ferromagnetic regime for strips with cylindrical boundary conditions. We 
expect a finite-size behavior of these three pseudo-critical temperatures of the type 

00 . 

v Q (q;L)=v c (q) + J2j^ (7-12) 

k=l 

where Q = E,C, U, and the {Aq^, ^k] are correction-to-scaling amplitudes and expo- 
nents, respectively. These quantities depend in general on q. We can try to estimate 
the critical temperature v c (q) by fitting the data to the above Ansatz (|7.12|) with only 
the leading correction-to-scaling term included 

v Q {q-L)=v c {q) + A Q L-^ (7.13) 

Although the exact phase transition temperature for the g-state Potts ferromagnet 
is known, it is useful to compare these pseudo-critical temperature with the exactly 
known values for the infinite lattice, since we do not know, a priori, which estimate vq 
will give the most accurate results. It is important to perform a comprehensive check 
of the method before studying a phase transition whose critical temperature is not 
known (for instance, the 3-state triangular-lattice Potts antiferromagnet; see below). 
We have done this and have found that the estimates coming from the pseudo-critical 
temperature ve are better by far than the other two; and the estimate v\j is more 
accurate than vc- More precisely, the difference between the extrapolated value for 
v c and the exact value is ~ 10~ 7 for ve for all the 1.5 < q < 4 values considered. 
However, for vy the discrepancy is of order ~ 10 -5 for q > 1.5 and ~ 10 -2 to 10~ 4 for 
q < 1. Finally, for vc the discrepancy is ~ 1CT 4 for q > 2, and ~ 1CT 3 for 1 < q < 1.5. 
In conclusion, we can establish the position of the ferromagnetic critical temperature 
v c (q) (over the whole range of values of q) by using the pseudo-critical temperature 
ve(L,L'). The results for vq and vu are at least two orders of magnitudes worse, 
and the accuracy also depends on the value of q: it worsens for q < 2 for vc, and for 
q < 1.5 for v v . 

One can try to extend the previous analysis to the antiferromagnetic regime. In 
Table |21 we show the estimates for the critical temperature of the 3-state Potts an- 
tiferromagnet using strips with cylindrical boundary conditions and widths that are 
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multiples of 3. This constraint is due to the (mod 3)-oscillations that appear in antifer- 
romagnets: in Figure ITBI we clearly observe such oscillations in the antiferromagnetic 
regime. Thus, we keep only the data with L = 3, 6, 9, 12 that is expected to be closer 
to the thermodynamic limit. 1 (When L is not a multiple of 3, the corresponding 
triangular-lattice strip with cylindrical boundary conditions is not tripartite, unlike 
those strips with L a multiple of 3 or the infinite triangular lattice; recall our earlier 
discussion of the chromatic number for these strips.) The value quoted on the row 
labelled "MC" comes from the Monte-Carlo study by Adler et al. [7S]. For L = 12 
we just list the estimate from the energy crossing since this is superior to the other 
two estimates, and this gives the value 

Vc (q = 3) = -0.796927(20) (7.14) 

The error bar quoted in ()7.14|) was roughly estimated by comparing the above result 
to the value of v c (q = 3) obtained by fitting the data points ve(L, L') with L = 3, 6, 9 
and V = L + 3, namely, v c (q = 3) = —0.796907. This is indeed a very conservative 
estimate for this error bar. Our results in (|7.14j) is in agreement with, and more 
accurate than, the estimate from [7K] listed above in eq. (|6.12|) . 
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lr ro compute the estimates ve(<] = 3, L,L') we need the transfer matrices T(Lp,g = 3). For 
Lp < 9, we have used the symbolic transfer matrices of section 4 evaluated at q = 3. For larger 
widths (i.e., L = 12, 15), we have computed numerically the corresponding transfer matrices at the 
particular value q = 3. 
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42 


132 


10 
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132 


429 


28 


24 
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429 


1430 


63 


49 
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1430 


4862 


190 


130 


9 


4862 


16796 


546 


336 


10 


16796 


58786 


1708 


980 


11 


58786 


208012 


5346 


2904 


12 


208012 


742900 


17428 


9176 



Table 1: Dimensions of the transfer matrix for triangular-lattice strips. For each 
strip width L we give the dimension of the transfer matrix for free boundary con- 
ditions N z ,tri,FF,L (which is equal to the Catalan number Cl), the dimension of the 
full transfer matrix for cylindrical boundary conditions (which is C^+i), the dimen- 
sion for cylindrical boundary conditions when translational symmetry is taken into 
account Nz,tri,PF,L, and the dimension when we project onto the subspace of reflection- 
invariant connectivities Nz, sq ,PF,L- 
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-0.7984897326 
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-0.7916555922 


-0.7942230532 
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-0.7971540641 




9 


-0.7952008646 


-0.7959971452 


12 


-0.7969905288 




12 






15 


-0.7969527708 




oo 


-0.79660 


-0.79668 


-0.796927(20) 




MC 


-0.79691(3) 


-0.79691(3) 


-0.79691(3) 



Table 2: Pseudo-critical temperatures for the 3-state Potts model in the antiferro- 
magnetic regime. For each value of the strip width L = 3k, we show the pseudo- 
critical temperatures vc, vu, and ve computed on strips with cylindrical boundary 
conditions. The row labelled "MC" shows the Monte Carlo estimate for v c (q). 
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Figure 1: Limiting curves forming the singular locus B q for the Potts model free 
energy for (a) v = — 1, (b) v = —3/4, (c) v = —1/2, and (d) v = —1/4 on strips with 
free boundary conditions and several widths L: 2 (black), 3 (red), 4 (green), and 5 
(blue). We also show the partition-function zeros for the strips L-p x (10L)p for the 
same values of L: 2 (□, black), 3 (o, red), 4 (A, green), and 5 (0, blue). Where the 
results for 2 < L < 4 overlap with those in ^ Ej (see also ^Hl EU E0| for v = —1), 
they agree and are included here for comparison. 
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Figure 2: Limiting curves forming the singular locus B q for the Potts model free 
energy for (a) v = — 1, (b) v = —3/4, (c) v = —1/2, and (d) v = —1/4 on strips with 
cylindrical boundary conditions and several widths L: 2 (black), 3 (red), 4 (green), 
and 5 (blue). We also show the partition-function zeros for the strips Lp x (10L)p for 
the same values of L. The symbols are as in Figure^ Where the results for 2 < L < 4 
overlap with those in |HE| and for v = — 1 (see also (SHIEIIBOJ for 3 < L < 5), they 
agree and are included here for comparison. 
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Zeros tri lattice q = 0, L = 2 p Zeros tri lattice q = 0, L = 3 p 
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Figure 3: Limiting curves forming the singular locus B, in the v plane, for the free 
energy, defined with the order f qn , of the Potts model for q = on the L F x oof 
triangular-lattice strips with (a) L — 2, (b) L — 3, (c) L — 4, and (d) L = 5. We also 
show the zeros of Z(G, 0,v)/q corresponding to the strips Lp x (10L)p for each value 
of L. 
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Figure 4: Limiting curves forming the singular locus B, in the v plane, for the free 
energy of the Potts model for q = 0.999 on the L-p x oop triangular-lattice strips with 
(a) L — 2, (b) L — 3, (c) L — 4, and (d) L = 5. This are essentially equivalent 
to the limiting curves for f qn at q — 1. We also show the partition-function zeros 
corresponding to the strips L-p x (10L) F for each value of L. 
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Zeros tri lattice q = 2, L = 2 F Zeros tri lattice q = 2, L = 3 F 
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Figure 5: Limiting curves forming the singular locus B, in the v plane, for the free 
energy of the Potts model for q = 2 on the Lp x oof triangular-lattice strips with (a) 
L = 2, (b) L = 3, (c) L = 4, and (d) L = 5. We also show the partition-function 
zeros corresponding to the strips Lp x (10L)p for each value of L. Where the results 
for 2 < L < 4 overlap with those in |1J Ej, they agree and are included here for 
comparison. 
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Figure 6: Limiting curves forming the singular locus B, in the v plane, for the free 
energy of the Potts model for q = 3 on the Lp x oof triangular-lattice strips with (a) 
L = 2, (b) L = 3, (c) L = 4, and (d) L = 5. We also show the partition-function 
zeros corresponding to the strips Lp x (10L)p for each value of L. Where the results 
for 2 < L < 4 overlap with those in |1J Ej, they agree and are included here for 
comparison. 
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Zeros tri lattice q = 4, L = 2 p Zeros tri lattice q = 4, L = 3 p 
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Figure 7: Limiting curves forming the singular locus B, in the v plane, for the free 
energy of the Potts model for q = 4 on the Lp x oof triangular-lattice strips with (a) 
L — 2, (b) L — 3, (c) L — 4, and (d) L = 5. Where the results for 2 < L < 4 overlap 
with those in jUE], they agree and are included here for comparison. We also show 
the partition- function zeros corresponding to the strips L-p x (10L)p for each value of 
L. 
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Figure 8: Limiting curves forming the singular locus B, in the v plane, for the free 
energy of the Potts model, defined with the order f qn , for q = on the L P x oo F 
triangular-lattice strips with (a) L — 2, (b) L — 3, (c) L — 4, and (d) L = 5. We also 
show the zeros of Z(G, 0,v)/q corresponding to the strips Lp x (10L)p for each value 
of L. 
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Figure 9: Limiting curves forming the singular locus B, in the v plane, for the free 
energy of the Potts model for q = 0.999 on the Lp x oop triangular-lattice strips with 
(a) L — 2, (b) L — 3, (c) L — 4, and (d) L — 5. These are essentially equivalent 
to the limiting curves for f qn at q — 1. We also show the partition-function zeros 
corresponding to the strips L P x (10L) F for each value of L. 
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Figure 10: Limiting curves forming the singular locus B, in the v plane, for the free 
energy of the Potts model for q = 2 on the Lp x oop triangular-lattice strips with (a) 
L — 2, (b) L — 3, (c) L — 4, and (d) L = 5. We also show the partition-function 
zeros corresponding to the strips L P x (10L) F for each value of L. 
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Figure 11: Limiting curves forming the singular locus B, in the v plane, for the free 
energy of the Potts model for q = 3 on the Lp x oop triangular-lattice strips with (a) 
L — 2, (b) L — 3, (c) L — 4, and (d) L = 5. We also show the partition-function 
zeros corresponding to the strips L P x (10L) F for each value of L. 
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Figure 12: Limiting curves forming the singular locus B, in the v plane, for the free 
energy of the Potts model for q = 4 on the Lp x oop triangular-lattice strips with (a) 
L — 2, (b) L — 3, (c) L — 4, and (d) L = 5. We also show the partition-function 
zeros corresponding to the strips L P x (10L) F for each value of L. 



58 



Energy tri lattice P bc's q=3 C H tri lattice P bc's q=3 




v 



(c) 

Figure 13: Thermodynamic observables for the 3-state Potts model on triangular- 
lattice strips of sizes L-p x oof- We show (a) the energy density E, (b) the specific heat 
Ch, and (c) the Binder cumulant U4 as a function of the temperature-like parameter 
v for several strip widths L: 2 (black), 3 (red), 4 (green), 5 (blue), and 6 (pink). — 1 < 
v < corresponds to the antiferromagnetic regime, while v > to the ferromagnetic 
one. 
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